NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



(NAS'l-C:;-1fi394 3) ^>>^A‘>AT. VOLUJIF 4: 

ATTITUDE DETBIMINAXIOJI ’liport (Jcit 

Propulsion L ib.) p Ui‘ A03/H1* A 01 

CSCL 0 ') B J iOi 1 > 


Seasat Final Report 

Volume IV: Attitude Determination 

Alfred J. Treder 


i 

V 

f 

V 

V 


July 1, 1980 



National Aeronautics and 
Space Administration 

Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 


— IT" iri'T-Tffi i *11 




JPL PUBLICATION 80-38, VOLUME IV 


Seasat Final Report 

Volume IV: Attitude Determination 

Alfred J. Treder 


July 1, 1980 


National Aeronautics and 
Space Administration 

Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 


The research described in this publication was carried out by the Jet Propulsion 
Laboratory, California Institute of Teohnology, under NASA Contract No. NAS7-1 00. 


PREFACE 


The Seas'it satellluu was launched at 01; 12: 44 GMT on 27 June 1978 from the 
Western Test Range at Vandenberg Air Force Base, Lompoc, California. The space- 
craft was injected into Earth orbit to demonstrate techniques for global monitor- 
ing of the dynamics the air-sea interface and to explore operational applica- 
tions. To achieve these objectives, a payload of sensors emphasizing all-weather, 
active and passive microwave capabilities was carried on the satellite. The 
mission was prematurely terminated on 10 October 1978 after 106 days of operation 
by a catastrophic failure in the satellite power subsystem. 

Major mission accomplishments were: 

(1) Demonstration of the orbital techniques required to support the 
mission and sensor operations. 

(2) Demonstration of the. simultaneous operation of all sensors for 
periods of time significant to global monitoring. 

(8) The collection of an important data set for sensor evaluation and 
scientific use. 

The early mission termination precluded: 

(1) Demonstration of the planned operational features of the end-to-end 
data system. 

(2) Collection of a global data set to meet overall geodetic and 
seasonal objectives and plans. 

This report, in four volumes, includes results of the sensor evaluations 
and some preliminary scientific results from the initial experiment team activi- 
ties. Scientific and applications studies will continue through FY 80, and will 
be included in the final version of this report. 


ABSTRACT 


On 7 July 1978 all sensors had been turned on, and Seasat began its 
successful mission. The attitude record from that date to the end of the mission 
has been archived on magnetic tape. In addition, a plot has been made of the 
available attitude data from each of the 1360 orbital revolutions, one plot per 
revolution, provided as a supplement to this report. 

Attitude behavior tended to be quite repetitive, revolution to revolution 
aiid day to day, for many days at a time. There were four major mission segments 
in this respect, each of which has been characterized by a single whole-revolution 
plot in pitch, roll and yaw. Each of these graphs shows the mean function of 
orbital phase, mean ± standard deviation of attitude variation, and the envelope 
of extreme values. The overall characteristics of attitude behavior can be 
inferred from these curves without reference to the individual revolution plots. 

The Seasat sun sensors constituted the only source of yaw attitude data. The 
sensor f ields-of-view limited direct yaw measurement to about half of the mission 
time, though each 6040-s revolution had at least a 900-s period of such measure- 
ment. Yaw data gaps caused by those limits were filled by use of a suboptimal 
interpolation algorithm. 

Of the four sun sensor heads, three observed the sun during the abbreviated 
mission and two of those showed misalignments on the order of 1/4 deg. Compara- 
tive analysis of data from one of the two misaligned heads produced an alignment 
calibration of accuracy sufficient to correct all the data from this head. The 
other misaligned head affected yaw accuracy during part of each revolution in the 
last nine days of the mission. 

The 3o accuracy of attitude determination (AD) was about 0.25 deg in pitch, 
0.16 deg in roll, and between 0.05 and 0.31 deg in yaw. Pitch and roll AD 
accuracy varied 25 to 50 percent with latitude, due to "cold cloud" error effects 
largely confined to lower latitudes. Yaw AD accuracy varied with sun position in 
the sensor f ield-of-view, cross-coupling of pitch and roll AD error into yaw, and 
alignment accuracy of the sensor head. Accuracy of yaw data gap interpolation was 
0.6 to 2,0 deg (3o) in the middle part of the mission, when the right Scanwheel 
experienced periodic sun interference, and the accuracy was 0.95 deg (3o) in the 
latter, quiescent part of the mission. 
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SIXM’ION 1 


INTIUmrCTION 


Tlus volume doeumonts ilu* aolual implementalion of iIk* Seasat Attitude 
Determination (AD) Hvsleiu and the eontents of tlu* attitude data files ^>;iaiera t ed 
by that system. Also included are the deviations from plan (*aused by tlie 
anomalous sun interference willi liorir.on sensors, infliplit c'alihration of sun 
sensor head 2 alip.nment and lun*izon sc*nsi^r biases, estimation of yaw interpola- 
tion [)arameters, sun and horizon sensor error sourci’S, and yaw inter])olat ion 
accuracy. lb\amples are p,iven of fliy,ht attitude data from all models ol‘ the 
Orbital Attitude Control System (OACS) , of the p,round proct*ss iiip, el fleets on 
attitude data, and td’ "cold cloud" effects on pitch and roll data. Altliouy.h 
not formally published, a supplemental set of plots, one for c>ach orbital revolu- 
tit>n, has been assembled for all definitive attitude data from 7 duly 1978 (day 
18S) throiudi tile end of the mission on 10 Octiiber 1978 (day 28d).* 

d'lu' Seasat attitude control subsystem controlled the. saLellilt^ ac’tively in 
pitch and roll usiny, ernu* siynals referenced to the horizon (see Section II of 
Volume 11). Yaw was maintained in a stable condition by ]) i t ch' momentum Mas that 
cross-coup 1 ed yaw momentum into roll momentum as the satellite movc'd around the 
orb i t . 


Knov\ 7 ledye of satellite attitude v/as required to a greater accuracy than 
control to properly locate the sensor measurements on the Earth *s surface and to 
provide Inputs for SASS data processing. The mission specification required; 

0,5 (dv\g (do) for control about each of the spacecraft control axes, and 0.2 deg 
(3:) for determination abcjut any axis. 

because of resource limits, the reference source selected for yaw was a set 
of sun sensors, which provided only partial coverage during the mission. The 
expectation was that the requirements could be met for all tliose times at which 
sun position information was available from these sensors. At other times, yaw 
attitude would be estimated by interpolation techniques, and the requirements for 
attitude knowledge would not be met completely. 

During the mission, several anomalies were encountered (sic Section VII of 
Volume II) that further affected the attitude control and knowledge accuracy. 

Because of these events and the importance of attitude knowledge to the user 
of the sensor data, a special effort has been made to describe the actual attitude 
history for the mission and to estimate the associated errors. The results of this 
work are. summarized in the following sections, and attitude histories are being 
made available as part of the basic data record. 


*In custody of Dr. Hiroshi Ohtakay, Guidance and Control Section (3A3), Jet Pro- 
pulsion Laboratory, 4800 Oak Grove Drive, Pasadena, CA 91103. 
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Ocher activities of this project are documented in separate volumes of this 
series ; 

Volume 1 Program Summary 

Volume II Plight Systems 

Volume III Ground Systems 

Abbreviations and acronyms used in this volume are defined in the appendix. 
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SKI' r I ON 11 


OKSOKimON OK An SYSTKM 


Thv SiMsat An svsUmu vKt'l H, lor juiriursi s o\ (his n»porl, 

oon.sislod ol that sali'llilc' hardw.irv' anO v;romul prvu'oioi soitwaio wluoh \no 
vUloU satollilo orionla*vIi>n inlormalion in a p.oooonlrio orbital nMiMouoo iraiuo. 
rho ininit.M to tlu* svstiMU woro llu* {uiu> Karth hoi'i;um, tuul sati'lliio onvironmoul, 
rho output vvMij a lU't oi t imo ordorinl l'^'limatos and rorri‘.spond inj\ unoorl a im i os 
ol Iho otlso' anp.lo:; Indwi'on tlu' Satollilo Allj»nmont RtMt'rom*o Axos (SARA') and 
tho p.oiu'ontrio oi'bital tolori'uoo axos. riu» svstiMu hardwaro oonsIstiHl ol at t i 
I uvlo siMisors u wo iulrarod hot iron juMnnoi\s tind two fU'l.s o\ .tfipoot sonsot s) 

and i ho i r oloolronios, ti'U'niOtrv livslom, and soioitoo inst rumont and tho .salt'h 
liti^ siruotiuv* on wluoh tho haixiwarv' wa?; mouttlod. Tho ^iot twaia' oon 

si.stod ol tho ni'liuitivo Attitudo Dot I'nuinat iiM\ tOAn) fuH Iwaro that providod a 
t imo historv ol tho attitudo o\ i ho .salt'll it roltuorn't' axoj; in t l\i' lotm oi 
dij\ital oomputor t apo lili'S* l.MSt' ]Uovidi‘d t lu' hardwaro aitd CVA'C doslj\nod, l>uill 
oporalod, .md maintainod t hi' pj'ound sot iwaro, Tho nAU t apo lilos woro input ti^ 
tho dlM. lootprint viol orminal (on soitwaro in suppvul o\ sc'ionlilio sonsvU' data 
proooss inp,. 


A. tTHh^niNAl'K SYSd’KMS 

l'ij»iu*o d I tRoliMS'noo ? ?) shows tho on orbit ooni i p,urat ion o\ tho Soas»u 
.satollilo. NvMuinal .*>paooorall axo;? i^Xoi Yo, anil Xo'l aro idonliliod in tlu* I 
nil' rolat ivo 1 tho 1 1 iphl path and Karth CnadiiA dirootions. 


1. Sati'llito Alipnmont Roloronoo Axos 

basoil vMi kiuuvModpi' al mas.s pri^port ios rolat ivo to tho Ap,ona i^oomotriv* 

.ixo.s uhrmipdt .spaoov'rall ilrawinps'), tho oriontation oT tho ovoratl .spaooor.ill 
l-rinoipal axos (pi'odiotovn was known to somo aoouraov. SARA (Ri'loronoo d \) ^ 

.1 vsuivoniont v^ptioal t v>o 1 roloronoo Iramo, was dot inod within tho .i I i j\nniv'nt li'ol. 
It w.o> dolinod bv a rot at ion I'o lat i vo to tho t ran.sl oi'rod .spav'ooral I roloi'onoo 
.ixos that alipnod SARA to tho prodiotod oriontation ol tho sisu'ooral l prinoipal 
.ixos. rho prodiotod prinoipal axt'.s wt'ix' t'xpoolod to 1 io wilt in O.d di'j\ ( VO o\ 
iho t ruo spaoovUMl t priiuupal axos. In any oaso, SARA wa.s tlu rtMorv'nm* Iramo 
t v> whiv'h all vUlu'r roloronoi' 1 ramos omboddod In spooilii' hardWtUt' woro laua- 
tion.illv alipau'vK In flipht, SARA was aotually iiuU't orminal o, but it ovmUi bo 
ii'v'roal »'vl with sullioii'nl aoouraov bv uso ol tlu' pround-^moasurod alipnnn'nt tram;* 
! v'rmal ivu\s. In tho oaso ol s\m sonsor hoad d, tins was aupmonli'vi bv Inllipht 
oalibration, Vvdiiv'h will bo disoussod in subiU'quont soot ions. 


ih*Mt Roloronoo Axt'S (tJoooonl r io , tlovulolio') 


Aii 

vlopoiuU'ni 


oharao t or I rod in Rol oronoo d- A, 

. t'.oiu'ral 1 V, tho X*-axis poIiUod 


t lu'so oov>rdinato .sv.'iliMiis woro vn’bil*’- 
In tho diroot ion v^l l lipht ami wa.'; 



Fi>\uro 


cumtalnod in tho orbit piano, tho Y-axLs pointed in tho dirootion opposite of 
the orbit momentum vector, and the Z-axis pointed in the direction of the Earth, 
For the geocentric orbit reference system, tlu^ Z-axis pointed to the center of 
the Earth (where the Earth rotation axis intersects the ecjuatorial plane), or 
the Z-axis was the negative of the orbit radius (position) unit vector. With 
the Y-axis parallel to the negative of the orbit normal, the X-axis was derived 
from the Y and Z axes, and was approximately in the inertial direction of flight. 
For the geodetic coordinate reierence frame, the Z-axis was pointing to the 
Earth, but it was normal to the local horizontal plane at the subsatellite 
point. The oblateness of the Earth caused a position-dependent angular offset 
between the geodetic and geocentric coordinates. 


B. ATTITUDE DETERMlNA'flON HARDWARE AND SOFTOARE 

Figures 2-2 and 2-2 (from Reference 2-2) shov; the locations and orienta- 
tion of the Scanwheel* horizon sensors and the sun sexisors on the Seasat space- 
craft. Tablt' 2-1 (from Reference 2-2) shovv^s the actual alignment of this 
hardware with rcvspect to SARA as measured on the griumd befcxre launch. These 
alignment values were taken into consideration in the ground processing of flight 
a L t i tiule t' e 1 erne t ry . 


1. Horizon Sensors (Scanwheels) 

References 2-1 and 2-2 give detailed descriptions of the design and nomi- 
nal operation of the Scanwheels, These devices provided the nadir reference for 
the OACS as well as the prime source of pitch and roll attitude data. 


a. Pitch and Roll AD Accuracy Augmentation . If the ht>rizon sensors 
were able to operate as a pair v^ith no deviation from nominal design parameters, 
they v^ould provide a geodetic reference to the OACS. Since the ground footprint 
software for science data processing requires attitude in geocentric coordinates, 
the least possible processing of pitch and roll attitude data would be the con- 
version from geodetic to geocentric coordinates. lUwever, pre-launch analyses 
of horizon sensor error sources (References 2-6 and 2-7) revealed the need to 
augment pitch and roll AD accuracy with additional ground processing. Accord- 
ingly, t'lSFC built several extra capabilities into their groxmd AD software, 
including horizon sensor bias determination (Reference 2-8), and correctitms 
(Reference 2-2) for biases and errors induced by systematic horizon radiance 
variations, orbital altitude variation, and Earth oblateness. 

Table 2-2 shows the results of the Scanwheel bias determinations performed 
during the mission. The first set of results was obtained from the mission 
period immediately after reliable OACS operation was established on the right 
(single) Scanwheel, during which time the prime Control Logic Assembly (CLA) 


^Tradename of Ithac.o, Inc. 
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Figure 2--3. roslllm\wS of Lho AtLltudo Dot enuinat ion and Control 
Hardware in the Spaeeeraft Body Viewed Along the 
Spaoeeraft Z-AxLs Toward the Karth 
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‘\\ nv*,l i' det in i t ^on^i : 

I'Uu'k ( t") a rv'tatiim about tbs C-tixis. 

Ceiu' t- > is a riUatlon about the rotated (due to t') Y-axis. 

IVist is a rotation about the rotated (due to x and iVaxis. 

The tvtations are -f* or - based on the rij^ht -hand- thumb rule where the tluimb is in the direction oT 
the vector and is a rotat ion in the direction of tlie fin^.ers. 

The rotations are from the nominal X,Y,Z spacecraft coordinate frame t»,; ,v semience). 

^'inflijdH calibrat ivm showi'd sip,nificant deviations from these values. 


I ' abU * 2 - 2 . Scanwhoel Bias Estimates 


Mission Period 
Kroni Which 
Rt'sultvS Weia’ 
Oht a ined 


Mission Period to 
Which Results 
Weri ‘ Applied‘S 


tU.A Pitwer 
Supply 
Sel ec t ed 


Fine Pitch Bias 

Kstimale, Uncort, 
dep, dep (TO 


Fine Roll Bias 

Kst imate, Uncert , 
dep dep ( 3 o) 


Pavs 18 h-.Mlh 


Days 1 H 8 -.M 1 


Days 


Days J 2 1-281 


0 . 12 


O.Oo 0.01 0.02 


0.11 O.OA 


0.10 0.01 


‘^FiU’ bias values used between days 220 and 22 b, see Section Ill-tO 
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pov/er supply 1 was operating eonf inuously . During days 220 and 221, sun 
interference with the Scanwheels resumed, and efforts to discover a favorable 
operating mode resulted in a permanent switch to CLA power supply 2. Hedeter- 
mination of biases subsequent to that switch showed an unexpectedly significant 
difference in n^ll, apparently caused by the switch in pov^er supplies. Defini- 
tive attitude processing corrected pitch and roll telemetry data for biases as 
a fvmction of the. CLA power supply, as shown in Table 2-2. Ri)th sets of biases 
were obtained from svm sensor head 1 telemetry data. 

Figure 2-4 shows typical factors that were added to pitch and roll telem- 
eti-y data to correct for the effects of Earth oblateness and systematic horizon 
radiance variations in both axes, ai\d for the effects of altitude variation in 
roll. This latter effect is due to orbital eccentricity and to differences 
between the mean orbital altitude and the altitude, implied by the bias voltage, 
which ill the GAGS replaced the left Scanwheel attitude signal. In the nominal 
dual-Scanwheel mode, such variations would have had no effect to first order, 
since the same effect occurred on both sensors and they would have been auto- 
matically nulled by the differencing performed in onboard roll signal processing. 
The single sensor mode in which Seasat was operated made roll subject to signif- 
icant altitude effects, as can be seen from the typical function (taken from a 
single revolution on 22 August 1978) shown in Figure 2-4. Actual correction for 
this effect was computed from the difference between 7168 Ion (3864 nm) and the 
actual satellite orbital radius; the actual radius was obtained from the defini- 
tive orbital ephemeris determined by (!SFC for eacb revolution. Therefore, the 
correction applied to roll for this effect varied in phase and amplitude from 
that shown in Figure. 2-4. Oblateness effects were also computed from definitive 
orbit data; these effects depended on the latitude of the subsatellite point, 
which is essentially normalized in Figure 2-4, so the actual oblateness correc- 
tion varied only slightly from that shown. 

Figure 2-.‘5 gives examples, for 3 months, of the corrections added to pitch 
and roll telemetry data for AD errors induced by variations in horizon radiance. 
These corrections are specifically for right (single) Scanwheel operation, and 
are significantly different in roll than would be necessary (and were planned) 
for nominal dual-Scanwheel operation. The analyses reported in Reference 2-2 
were reworked by GSFC immediately after the Seasat launch, vdien it became appar- 
ent that the nominal mode could not be used. This report, and the GSFC mission 
report, are the first documentations of the results of those reworked analyses. 
All definitive attitude files released in final form by GSFC were corrected with 
factors appropriate for single Scanwheel operation. 

b. Example of Pitch and Roll Flight Data . Figure 2-6 shows, as an 
example, pitch and roll data from a single revolution (172) from day 190 (9 July 
1978), both with and without corrections. Both plots were made from data at 5-s 
intervals, but a 5-s running average filter was used on the uncorrected graph, 
and a similar 120-s filter on the fully corrected plot. The latter mode was 
used for all definitive attitude files from day 188 through the end of the 
mission . 
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Figure 2-5. Horizon Radiance Corrections 



Corrected and Uncorrected Examples of Pitch and PvOl 
Flight Data 




Tht' of pitch hiius and ohlatcncos corrections arc easily seen from 

compai isons of Fi^nuvs 2-^ and ll~6; this comparison for roll is not as obvious 
since the correction for altitude in lv4);ure is not the correct phase fm* 
rev 172, but the oblateness effect is fairly cliMr. "th'tld i'loiui"' effects are 
discussed in detail in Refiu'cnce 2-9 but, as an example, note the positive 
0*08-deg trian^uilar function that occurs simultaneously in both axes at about 
6100 s of day; this ivS the effect of a small cold feature v>bsei'ved by the 
threshold-setting", portion of the trailing edge field of view (FOV) of the right 
Scanwheel. The 1 20-s nirm'ag average has reduced the rcvSnlting. iittitude deter- 
mina l i i m e r r o r 1 r om 0.08 d c g t* o 0.0 1 d e g , a re du c t i on o f ab ou t 60 e r c en t . 
Shorttu* duration features would be filtered out more heavily, v/hile longer dura- 
tion features v;ould have imly their high-frequency components smoothed out. The 
effects on AD accuracy of the ninning average filter are discussed later. 


2, Sun Sensors. 

The Adcole sun sensor system consisted of two sets, each set containing two 
2-axis sun sensor heads (Figure 2-7) and an electronics package to condition the 
sensor signals for digital telemetry. The resolution for each sensor axis was 
1/256 deg, and the absolute accuracy was 0.05 deg (3o) relative to the alignment 
mirrors. These sensors were not a part of the GAGS, Ciid were used only for yaw 
at-titxide determinat ion and horizon sensor bias determination. 


Sim Sensor Coordinate System. The sun sensor coordinate system 
(Reference 2-2) was defined so that the positive Z^-axis was directed along the 
optical boreslght of the sensor. The and axes defined a plane normal to 
the Zj^-axis and were oriented as follows. The transformation from the SARA to 
the sensor axes was defined by an ordered 3-2-3 Euler rotation: clock angle 

rotation of the sensor about the SARA positive Z-axis, followed by a cone angle 
rotation about the new positive Y-axis, followed by a tx^ist angle of -90 deg 
about the new Z-axis as shoim in Figure 2-8. The difference between the measured 
twist angle and -90 g is defined as the rotation angle. The Xcj^-axis is 
coll inear with slit v v>f the two-axis sun sensor, and the Y^,-axis is col linear 
vxrlth slit A of the sun sensor iis slunm in Figure 2-9. Pre-launch measured align- 
ment angles of these sensor axes x^7ith respect to SARA are given in Table 2-1. 


8 tm Sensor Field s -of-View . Figure 2-10 shows the FOV for each of 
the four sun sensor heads In the system. The coordinate system is spacecraft 
SARA, in a polar presentation with xenith (180-deg cone) at the center. The 
locus of nominal sun positions is given at intervals of 4 deg in sun beta angle 
over the ctmrse of the mission. Note that sensor head 4 never saw the sun, 
because the mission ended prematurely. The locations of the heads were chosen 
during the design phase to; 

(1) Maximize the mission time for which the sun would be visible 
to the sensors , 

(2) Ensure that every orbital revolution could have some sun 
observab i 1 ity . 
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ALIGNMENT MIRROR REFERENCE 


Figure 2-7. Fine Sun Sensor Head (Containing Independent 
One-Axis Sensors) 
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(3) Ensure a small FOV overlap between heads, where practicable^ 
for relative calibration purposes, 

(4) Keep all spacecraft structure out of the 64- by 64-deg FOV to 
avoid any interference or glint problems. 

The final selected orientation met all of these criteria in a suboptlmal but 
acceptable manner. 

Figure 2-11 presents nominal sun sensor coverage as a function of mission 
time. Coverage is defined here as the time from ascending node and the subsatel- 
lite latitude over which a sun sensor (Identified by number) could observe the 
sun, given the nominal 64- by 64-deg FOV, no significant attitude deviations, and 
expected sun beta angle at the time. Note the slight overlap of brads 1 and 2 
near day 240. The meandering double line curves denote the times when the sun 
was within 1.5 deg of the right Scanwheel FOV, which will be discussed later. 

In flight, the sun sensor FOV had to be cut off at 30 deg from boresight, 
leaving a 60- by 60-deg FOV. This was done by setting appropriate limit para- 
meters in GSFC ground software. The cause for this action, which significantly 
reduced the mission time of sun observability, was that the sun sensor system 
produced telemetry signals for sun positions outside the ±32-deg boundary, and 
these signals were indistinguishable per se from honest readings for sun posi- 
tions within the 30- to 32-deg band. That this could happtni was apparently 
understood by LMSC and Adcole early in the design phase, but x^ 7 as not adequately 
communicated by them to JPL and GSFC, As a result, GSFC software as delivered 
incorporated no procedure for resolving the data ambiguity. Wien the problem 
was discovered, no time or resources could be committed to make the required 
software changes, so amputation was accepted as an expedient fix. It was still 
possible to reclaim this data through expert intervention by the software 
operator, but this resource-intensive method was used only for special occasions, 
such as for saving the FOV overlap region betweens heads 1 and 2 for calibration 
purposes. Figure 2-10 shows both the nominal 64- by 64-deg FOV and the 60- by 
60-deg FOV realized in flight. (See Section 3,1.1 of Reference 2-9 for a 
detailed discussion of this problem.) 

Actual sun coverage was affected by attitude deviations near the FOV 
boundary and by loss of telemetry data for any number of reasons. Figure 2-12 
shows the actual sun coverage boundaries experienced on the mission, plus many, 
but not all, data gaps due to other causes and an outline of predicted sun sensor 
coverage. Figure 2-12 is provided to help identify the size and location of all 
major yaw attitude data gaps. Figures 2-11 and 2-12 can be compared by noting 
that the first quarter of each revolution is shown at the top of Figure 2-1] and 
at the bottom of Figure 2-12. Except for scale, all other factors are the same. 


c. Yaw Determination from Sun Sensor Data . The algorithm used by GSFC 

to determine yaw attitude from sun sensor data is described in Reference 2-2. 
Digital sun sensor telemetry is converted to cone and cross-cone aspect angles 
(see Figure 2-9) by applying a polynomial function derived from ground calibra- 
tions. These two angles define the actual sun vector in sun sensor coordinates. 
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Figure 2~11. Nominal Sun Sensor Coverage 
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The diffurt'nooK betwoen tbi.>. actual sun vector and the vector expected at that 
orbital position are assumed to be caused by attitude perturbations. Using 
pitch and roll data estimated from Scanwheel telemetry permitted a closed-form 
solution for yaw. The accuracy of yaw determination was, therefore, dependent 
on pitch and roll determination accuracy, and the positiini at that time of the 
sun within the sun sensor FOV. 

l-Hien no sun sensor could see the sun, yaw could not be determined directly. 
An indirect method of interpolating these yaw gaps was used if telemetry data 
(including roll attitude) were available. This method will be described in the 
fol lowing paragraphs. 


3. Yaw Data Gap Interpolation 

The following information is adapted from Reference 2-10. 


a. S ummary of Method . Yaw data gap interpolation involves the problem of 
filling out the record of a real function of time, the measurement of which is 
desired to be continuous but which is arbitrarily discontinuous in its physical 
Implementation. A typical application is completing the attitude data record of 
a nadir-pointing spacecraft in low planetary orbit where a sun sensor is chosen 
to be the attitude reference, in one axis. Such a system will have attitude data 
gaps due to sensor FOV limits and occultations of the sun by the spacecraft and 
the planet. In the case of Seasat, maximum science return from the mission 
required complete attitude knowledge at every point in the orbit to an accuracy 
level 2.5 to 4 times better than tne control of attitude. To produce continuous 
and accurate attitude determination by addition or substitution of star sensors 
would have been too expensive, as would have been the development and continuous 
use of an adequate spacecraft simulation on ground computers. Instead, a more 
limited amount of mission funds was devoted to the support of attitude determina- 
tion, and those funds were used to extract as much attitude information (and, 
therefore, scientific return) as could be obtained under the circumstances. 

This is a typical design trade-off, leading to a suboptimal but suffi- 
cient attitude determination system. The method developed to provide a reasonably 
good solution to the constrained data gap problem is generally applicable to any 
system characterlaabl e by a Gauss-Markov process, and it produces a stable, 
closed-form solution under all conditions. 

The general problem can be stated as follows: given the data record 

D(t) of a system parameter Y where at arbitrary times within the range, of this 
record no direct observations of Y exist, provide sufficiently accurate estimates 
of Y using all available types of Information that pertain to Y and describe the 
variance of the estimates. 


The desired solution in the widest sense would provide estimates 
with variance less than or eqvial to the variance of directly measured data. 
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The practical solution given here will provide estimates at all request times, 
the variance of which is maintained as close to 0 ^ as the (variable) quality 
of input information will allow. The algorithm can be incorporated into a data 
processing system with little or no risk, and can be invoked reliably and 
automatically on a regular basis at a fixed, reasonable cost in throughput time. 
The products of the algorithm have predictable bounds of uncertainty. The 
algorithm was designed generally enough so that even after it was coded and 
incorporated into the data processing system, its output could be upgraded by 
making appropriate adjustment of its input parameters. The amount and quality 
of such algorithm support was dictated by cost constraints and the availability 
of the appropriate information. 

The method uses three classes of information: 

(1) Direct Indicators of the system parameter Y when they exist, 
the measurements of which comprise the data record D(t), and 
the statistical extension of this record into time intervals 
of no direct measurement (data gaps) using a modification of 
standard autocorrelation techniques. 

(2) An a priori predictor Yp of the behavior of Y, based on the 
most accurate available model of system behavior. At the 
least, Yp would describe only the design control limits of Y 
about its nominal value, and at the best would be a complete 
system model capable of exact prediction. 

(3) One or more indirect indicators of Y based on observations 

of other system parameters. These functions may be used even 
if their errors are partially correlated. 

The algorithm first extrapolates the available direct data D(t) into 
the data gap from both edges, de-weighting it exponentially as a function of 
distance from gap edge. The a priori predictor function Yp is given increased 
weight as the extrapolation is de-weighted. The variance of this combined 
estimate is calculated on the basis of the stochastic process (first order Gauss- 
Markov) assumed for the variation of Y with respect to Yp , the autocorrelation 
constant determined from the observed behavior of (Y-Yp) , and all of the appli- 
cable measurement, noise, and parameter variability factors estimated a priori , 
outside the algorithm. The indirect measurement of Y (class 3, above) is 
treated as an estimate with errors that may be correlated with those of the 
extrapolation estimate. The correlation is estimated a priori (from previously 
observed behavior) and then used in combining the estimates in the linear mean- 
squares sense. The variance of the algorithm result will lie between the 
variance of the linear mean-squares optimal combination of estimates, as the 
anticipated minimum, and the variance of the most accurate component of the 
combination at the local time of the estimate, as the maximum. 
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b. Interpolation Teclmldue.s . Data interpolation is an old problem, 
and many techniques have been evolved. The following is a list ol approaches 
that could be applied to the arbitrary data gap problem: 

(1) Data Smoothing . iNllien the time spanned by the periods of 
relatively continuous data is an appre«. .able, fraction of 
the system time constant, and the periods of data outage 
are short in the. same frame of reference, one of the many 
smoothing techniques can be used. Fourier series approxima- 
tions and polynomial or spline fits to the data are approp- 
riate techniques under the limiting conditions. As the 
data gaps get larger than the periods of observations, 
these, techniques must be limited more strongly to their 
lower order terms, and accuracies drop rather rapidly. 

(2) Extrapolation by Autocorrelation . To the extent that system 
behavior can be modeled via a stochastic Gauss-Markov process, 
the value of the system parameter at times greater than the 
last observation (and at times less than the next observa- 
tion, by symmetry) can be inferred from the parameter's 
statistical correlation with itself. The true process 

must be nearly wide-sense stationary in the local time 
frame for meaningful estimates to be made, but the error 
caused by relatively small deviations from this condition 
can be handled in many cases by tracking the changes in 
autocorrelation and by appropriately scaling the estimation 
uncertainty. As time from the last observation increases, 
the information in the estimate decreases exponentially, 
and, therefore, in a practical sense its value disappears 
rather quickly. 

(3) A Priori Prediction . A function or set of functions of 
time is used to predict the behavior of Y, based on a 
characterization of the system design and a prediction of 
the disturbances. If resources permit, a complete set of 
system equations is modeled, taking account all system 
internal characteristics and responsive to all external 
drivers and disturbances. Such a model is often expensive 
to construct and maintain, but can usually be made as 
accurate as desired. Less complete models are capable of 
less accurate results, given that no observations are noise- 
free. The minimum such model is an estimator of the mean 
behavior of the system parameter over all time. Given no 
observations of the system in operation, the maximum error 
of this estimator is the control limit of the system in its 
operational environment. 
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(4) Indirout M easurement . If lliero uve othor sysuoni parameters 
tiiat are observable when the parameter of fnterest is not, 
some form of a system model may be used to relate the 
behavior of these parameters to the behavior of the desired 
function. A known eross-eoupl in^» of axes can be useful, as 
can a statistical correlation with the Y parameter that is 
believed strong in the local time frame. As an example, a 
nadir ^-pointing orbiter has an Inertial cross-coupl ing between 
the two spacecraft axes (roll and yaw) in the orbit plane. 

The accuracy of such an estimate might be much less than 
that of a direct measurement and yet be xxseful, since it 
contains information near the time of interest. It is, 
therefore, more likely than the previously listed methods 
to show the results of random disturbances occurring in the 
data gap. 

The hybrid method presented here to handle the arbitrary data 
gap problem includes aspects of all of the above techniques. 
When data gaps are short relative to the autocorrelation 
time (t) of Y, the method becomes basically an exponential 
process of high accuracy. As gap length approaches and 
exceeds i, the a priori function and indirect measurements 
have stronger effects on the total estimate, and eventually 
become predominant. In all cases, however, the estimates 
near the gap edges approach or exceed the accuracy of 
directly measured data, with uncertainty growing with 
distance from the last observation to error bounds that 
remain within system limits even for very large gaps. 

Unless the chosen system model requires it, no mode switches 
need be thro\m to effect this broad range of performance. 

At each time point within the gap, the quality of the esti- 
mate is a fair approximation of the input information 
quality at that time. 


c. Application to Seasat . The Seasat mission that began in vlune 1978 
used the algorithm described below in its ground support of attitude determina- 
tion. Seasat was a nadir-pointing Earth-orbiter with momentum bias perpe^ndi- 
cular to the orbit plane (pitch). Control with respect to local nadir was 
obtained with horizon scanners, which were also the primary sources of attitude 
knowledge in pitch and roll. Rotation about the nadir (yaw) was limited by the 
stiffness of the momentum bias and the inertial cross-coupling into the controlled 
roll axis. Roll was, therefore, an indirect measure of yaw under steady-state 
conditions, although actual conditions made roll a relatively poor yaw indicator. 
Roll rate was numerically estimated from roll position data, and it was sometimes 
a relatively good indirect measure of yaw disturbances through the gyroscopic 
effect of the momentum bias. Direct measurements of yaw were obtained from the 
sun sensors. The length of the yaw data g£ip varied from 0 to 80 percent of 
orbit period in a relatively predictable pattern (Figure 2-12). Full computer 
simulation of the relatively complex Seasat control system and its environment 
for AD purposes was ruled out by cost constraints. Observation of yaw during 
full-data periods of the mission permitted a limited-capability empirical yaw 
predictor, and evaluation of roll and roll rate as indirect measures of yaw. No 
assessment of autocorrelation and cross-correlation was attempted. 


d. 


Definitions . 


t 




t2 

Y(t) 

Yu(t) 

Yp(t) 


a 


2 

c 


^1 

■^2 

R(t) 

R(t) 


Time from ascending node (the point in the orbit where the 
subsatellite point crosses the planet equator in a northerly 
direction) . 

Time (in the above sense) at which the flow of directly 
observed data is interrupted; i.e., initial time for the 
data gap. 

Final time (in the above sense) for the data gap. 

True yaw behavior as a function of t. 

Directly observed yaw data, as a time sequence. 

Function of t which predicts a component of yaw behavior. 

For Seasat, the following algorithm was used: 


Yp(t) 


4 

= + 2 ] K^cos( 27 Tit/P^ + X^) 

i=l 


where are scaling and phasing constants chosen empir- 

ically to minimize o2, and Pq is the orbital period. 

Variance of Yp(t) with respect to Y(t). The value chosen 
must be consistent with observed yaw behavior, the effects 
of predicted disturbances not modeled by Yp(t), and 
predicted system- imposed limits on Y(t). 

Variance of observed yaw data v/ith respect to Y(t); i.e., 
absolute measurement accuracy. 

Autocorrelation time constant of pre-gap yaw data. 
Autocorrelation time constant of post-gap yaw data. 

Roll attitude data at time t. 

Roll rate numerically derived from smoothed roll attitude 
data near t. 
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Y^(t) Estimate of yaw at time t based on observations of roll 
and roll rate, computed throi4*;h a model of the yaw/roll 
eross-’coupling. The form of this set of effects wa^s 
assumed for this analysis to be: 

Y 3 (t) = Kyi^ R(t + tj^) + Kyrd R (t + 


where Ky^, KyRn ^caliiiK constants, and Lr, tRn are 
phasing’ constants, chosen empirically to minimize 

Variance of Y3(t) with respect to Y(t). It is evaluated as: 


^I3M* 


E 



(L) - Y.3(t) 



9 

less the measurement noise o^, where the expectation is 
taken as ensemble over multiples of orbital periods* 


Healing factor and rime constant for assumed model of 
cross-correlation of errors in and 'in the interpolation 
estimate Y^ (defined in text below) . 


To begin the construction of the estimation algorithm, treat the function 
Y(t) as a stochastic process or sequence. Assert for these purposes that the 
deviations of observed data Yj^Ct) with respect to the a priori deterministic 
function Yp(t) are random or of unknoTO causes. This models Y(t) as the sum of 
Yp(t) plus noise: 

Y = Yp + Yr 

where Yj- is a random variable modeled by a Gauss-Markov process having 
the following properties: 
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By the Gauss-Markov assumption: 


p(At) = 

The equation for the measurement of Y^(t + At) by the observation Y^^Ct) isi 

YhCt) - Yp(t) = p(At) . Y^(t + At) + n 

The expectations given above yield: 

2 2 ^ 2 . 2 

"n ■ % * > "c 

Then the estimate of Y(t + At) obtained from Yj^(t) is: 

Y(t + At) = Yp(t + At) + Y^(t + At) 


p [Yj^(t) - Y^(t)] 


2 . 2, 2 

Op + (1 - P ) 


= Y (t+At) + 7 

P ' 


-1 + _ 
2 2 

0 

L c 


4. , y - 2 V 2 

»i, + a - P ) 


2 2 2 

Y c 

r 


2, 2 

Op + (1-P ) O^ 


2.2 

0^ + 0 
D c 


where p in the above equations is the previously defined p(At). 

In the application of this estimation process to the data gap problem, 
there are two Yp observations, one at each end of the gap. We can take advantage 
of this extra information by expansion of the previous expressions. 


Let 


Pi = ^ 


-|t-tj /t^ 


and 


P2 = e 


-|t-t2|/T2 


t^ < t < t2 


( 1 ) 
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Then 






2 


D 


o-^ll 




2 

'2 


2 

c 

■^2 

(1-P2I 

2 

n 

C 


-1 


( 2 ) 


Y^(t) = Y (t) + a (t) 
I pi 




I' D' r 




2 j. /i 2 \ 2 

'n + (1-^ i) % 


’I) + (1--2) 


2 

c 


( 3 ) 


where the subscript I is used to indicate interpolation. 

In the limit as (t2 - t^) becomes much less than t j and Yj(t) approaches 
a straight-line interpolation between Yp(tj) and Yp(t2). In the middle of data 
gaps much larger than (tj + T2)> Yj(t) approaches the value of the a pr i ori 
function Yp(t)j and the variance of YpCt) relative to Y(t) approaches , which 

is the variance of Y (t) relative to Y(t). 

P 

The Yj(t) estimate reflects the Information contained in the data at the 
gap edges and in the a priori function Yp. The set of indirect measurements of 
Y(t) expressed by the previously defined function Y3(t) can be used to enhance 
the accuracy of Yj(t). The degree of enhancement is inversely proportional to 
the amount of statistical correlation between Y3 and Yj. 

A combination of estimators with correlated noise may be formed optimally 
with respect to their mean squared errors by using a set of multipliers that 
satisfy the orthogonality criterion. For the set of estimators Yj and Y3, the 
optimal combination estimator would require the characterixation of the variance 
of Y3 and the covariance of Yj with Y3 as a function of time. The estimator is; 

Y(t) =L^(t) Y^(t) +L2<t) ( 4 ) 

oj (t) = L^(t) Op(t) + L^Ct) a^Ct) 

o 

PfsCt) 0^(t) 


+ 2 L^(t) L^Ct) 


( 5 ) 


where 


Pl3(t) 


E I Y(f.) 


E LY(t) 


- Y,j(t)] [Y(t) - YgCt.)] } 

- Y^(t)rf. E I Y(t) - Y3(t)rf 


0 (t) - cf (t) 0 (t) 

^ _ 2 i d . 

^ ^ ^3(t) 


L 2 <t) = — 


^iCt) - Oj(t) o^Ct) 


o^(t) + cT^Ct) - 2 "3^^^ 


In many practical applications, and in particular for Seasat, a rigorous 
characterization of as a function of time is not necessary for suboptimal but 
adequate estimation accuracy. This is particularly true when the variances of 
the estimator components are widely disparate, in which cases the combination 
estimate strongly favors the more accurate component regardless of correlation 
values. For this reason, equal to zero for Seasat. 

2 

If Y 3 (t) has a sufficiently constant variance ( 03 ) relative to Y(t) over 
the local time domain of Interest, 03 may be modeled as the ensemble average 
over that domain. This simplification is an economy in the design of the com- 
bination estimator, and was implemented for Seasat. Because of the simplifying 
assumptions on 03 and Pi 3 (t), plus the possible modeling error in the Yj estimate, 
the variance calculated by equation (5) is expected to be an optimistic esti- 
mate of the true estiinator variance. The expected upper bound on estimator vari- 
ance is O 3 or o^, whichever is smaller at the time of the estimate. 

Equations (1) through (5), wita the associated definitions of parameters 
and qualifying statements, constitute the yaw interpolation algorithm for Seasat. 


e. Algorithm Performance on Simulated Data . Figure 2-13 shows some 
typical design simulatima results for Seasat steady-state attitude in the roll/ 
yaw plane. Roll, roll rate, and yaw data are displayed. From this data, other 
similar simulations, and knowledge of Seasat design, the parameters of the yaw 
estimation algorithm were estimated and are listed in Table 2-3. Because thes^ 
values were based on analysis of only a limited number of simulations and no 
flight data, they were only preliminary for Seasat purposes. However, they were 
appropriate values for demonstrating the operation of the yaw estimation algor- 
ithm and its effectiveness on available simulation data. 
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UvSod to I'hook out tho ooiu'o|U . Fipau'i' d-lA shows tho otUimalor oom]umouts with 
notation as dol’inod provio\uily lor l lu' alp.orithm applied to the data of FIp,*^ 
ui‘o x\n artilioial yaw data p.iip Wiis oroatod hotwoon tho time poinlt; marked 

tj and ti, and the Y data points made unavailable to the estimator are indieatiHl 
hy a d ashed line, NiU e that no attempt was made lo lit sK>pe and aeeel oration 
at the boundariet?; thoupji unaesthetie and contrary to system dynamics, the 
deviations causeil hy these breaks are properly accotmtinl in the Vviriance est !-• 
rated by the al^^uMthm, il' the model asstimpt iom; are correct , 

Pip, lire *''la also shows the optimistic bound of pretlictc‘d estimator error 
(solid i'urve, uj^per ripditl at \o and the actual deviation of the estimate from 
the thidvKnO datvi versus time lor the example. This is for illustration only, 
since a proper comparison would involve the ensemble ol estimator deviations 
over a larpe number of sample data paps. Such a comparisvni is made later. 

Fipui'o illustrates alporithm performaace on dat:i pa]is ranpinp irom 0.7 

lo 1.8 autocorrelat ion time Cil lenpths. Fipure I’-lb illustrates perlormance for 
data paps of about b; , a situation that faced the estimation alporithm about SO 
percent of the time it was used for Seasat (see Fipure d-lH. 

Overall performance of the majtu‘ yaw estimator compinients is piven in 
Fipure d-1 7 as three times imis error of each component with respect to true yaw, 
averaped over all diita points from simuhited paps of the same lenpth. The data 
paps were art if lea lly created in a xmifoimi manner as a prid on the available 
simulation data, as a quas i-*Mont i Ckirlo approach. All yaw behavior as represented 
by these particular simulations was, theivfore, piven approximately the same 
opportunity to affect these performance statistics. Note the relatively stronp 
pertornumce of the interpolator (iSyj) for shorter data paps, with error prowinp 
to apprcMch that of the priori function for larper paps, as would be 

expected. Roth the a pr iori prediction O’lO aad the cross-corre lat ion term (Y^)* 
have aiqu'v^x imat e ly constant statistics; these results are forced by tlu' lUMiimira- 
t ion of and perfiumied before these sample statistics \yfovo taken. Note, 
finally, tlxit in the bOO™- to bOOO-s pap replon the total yaw estimate statistics 
(38^) are always better than or equal to the host of its componmUs, and not as 
ptufd as the variance of the optimal combination estimate as predicted by the 
alporithm. The interestinp dlverpence from this condition fiu' shorter dO‘:x paps 
is due partly to the small number of samples that went into the statistics, and 
partly to the inaccuracy of lonp-*scale ensemble averapes of and e*; with 
respect lo these particular samples. 

The final comparisim is offered v>rith a note of caution. Tlu' rms error of 
each estimator ctmiponent was noniuilined to the rms value of yaw data itself for 
each sample data pap of the same lenpth, and the avempe of these rat ios is 
presented in Fipure 2--I8 as a fimct ion of data pap lenpth. These comparisons 
are measures of the improvement offered by components of the yaw estimation 
alpm'ithm over the rather simpler method of estimatinp yaw 1 be aero everywhere 
in the pap. It is an interestinp comparison, and shows the sipnificaiU improve- 
ments available, but the linearity of this relationship needs to be demonstrated 
on a wider ranpe of sipnal amplitudes. More comprehens i ve est imati-s oi perform- 
ance level on flipht data will be piven in subsequent sections. 
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YAW ANGLE, deg 



Figure 2-14. Typical Estimator Components 




YAW ANGLE WITH RESPECT TO DATA NULL, deg 



Figure 2-15. Typical Simulated Yaw Data Record, Example 1 
(Data Gaps Filled With Estimates) 
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YAW ANGLE V7/RESPECT TO DATA NULL, deg 


0.20 



Figure 2-16. Typical Simulated Yaw Data Record, Example 2 
(Data CJaps Filled With Estimates) 



ESTIMATION RMS ERROR, deg 
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Figure 2-17. Yaw Estimator Performance Accuracy (on Simulated Data) 
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DEFINTTIVK ATTITUDE FILE 


1. Introduction 

The Definitive Attitude File (DAF) is a record of the best estimates of 
spacecraft attitude defined on geocentric orbital reference axes as a function of 
time. The DAF was produced by Code 580 at GSFC and was recorded on an Attitude 
Orbit Tracking (AOT) tape generated for each day of the mission by the GSFC 
Information Processing Division (IPD) . That tape was delivered by the Instrument 
Data Processing System (IDPS) at JPL. The DAF software capability was operational 
before launch on 27 June 1978, and the IPD capability x?as operational somewhat 
later. The functional flow of attitude data at GSFC is sho^m in Figure 2-19. 


2. Functional Description 

A file of determined satellite attitude spans the same satellite data day 
as the contents of the Telemetry Master Data File (TMDF) . The file begins at 
00:00:00:000 and ends at 24:00:00:000 GMT of the data day, according to require- 
ments. Actual files deviate somewhat from this ideal, with the deviations always 
in the sense of less data. 

Each data point of determined satellite attitude is expressed as a time 
(milliseconds of day), and the set of Euler angle rotations which at that time 
will convert to the SARA any vector expressed in the geocentric orbital reference 
axes defined as: 

(1) Axis 1 (X): in the inertial direction of flight (derived from 

Y and Z) . 

(2) Axis 2 (Y) : parallel to the negative of the orbit normal. 

(3) Axis 3 (Z) : toward the Earth center of mass. 

The Euler angle rotations are in the 3, 1, 2 order of axes, which corresponds 
to satellite yaw, I'oll and pitch. 

Yaw attitude at times the sun sensor can see the sun and pitch and roll at 

all times (subject to telemetry availability) were to be determined with a total 

absolute accuracy of 0.17 deg (3o) for each axis, provided that all spacecraft 
error sources were meeting their attitude determination and control requirements. 
In flight, all attitude data was obtained in a non-nominal control mode, which 
degraded the accuracy of Scanwheel telemetry data. Therefore, although DAD soft- 
ware was designed to obtain the highest accuracy realizable from the data, the 
basic data quality was not good enough to permit attainment of the accuracy goal. 

Yaw attitude at times the sun sensor system could not see the sun was 

determined with the yaw estimation algorithm described in the preceding para- 
graphs. Yaw determination accuracy during those times is whatever is provided 
by the algorithms . 
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SEASAT 



Figure 2-19. GSFC Attitude Determination Data Functional Flow Diagram 






Tho. time and frequency of the attitude points were selected to be synchro- 
nous with the Definitive Orbit File (DOF) data points. That is, since DOFs were 
provided on even minute time marks, the DAF contained data points at each of those 
times, and other attitude points were spaced between those times at 5-s intervals. 

During periods for which no spacecraft telemetry exists in the TMDF, atti- 
tude points were not required, although they may have been included on the file. 
There were no zero-filled attitude logical records; i.e., data gaps were com- 
pressed out. Gaps are defined as time intervals when no telemetry exists. 


3. File Organization and Contents 

Each physical record in each DAF was 3476 bytes long. The first record was 
a header record. The remaining physical records each contained 144 attitude data 
sets. The organization of the header is shown in Figure 2-20, and its contents 
are listed in Table 2-4. The organization of the data records is shown in Fig- 
ure 2-21, and their contents are listed in Tables 2-5 and 2-6. 

Each physical record contains 144 logical records. There can be up to 121 
data records for a 24-h period. The last record in such a case contains only the 
one data point at 24:00:00, and the rest of the record is filler. 
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DATA BASE IDENTIFICATION 


OUTPUT REPETITION PERIOD (s) 


MAX PITCH (deg) 


MEAN PITCH (deg) 


MIN PITCH (deg) 


MAX ROLL 


MEAN ROLL 


MIN ROLL 


MAX YAW 


MEAN YAW 
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DAY OF YR 
MSEC OF DAY 


YAV/ FLAG 


E 3 . BY 


MSEC OF DAY 


PITCH FLAG 
^YR 


ROLL FLAG 


• LOGICAL RECORD 1 


DAY OF YR 
MSEC OF DAY 


YAW FLAG 


MSEC OF DAY 

PITCH FLAG 

ROLL FLAG 

^YR 

o-Y 


> LOGICAL RECORD 2 


DAY OF YR 

MSEC OF DAY 

MSEC OF DA” 

PITCH FUG 

ROLL FUG 

YAW FLAG '’YP 

'’YR 

cty 


LOGICAL RECORD 144 


FILL ZEROS (20 BYTES) 


Figure 2-21. Definitive Attitude File, Data Records 
(Physical Records 2 through N) 


0 


Table 


Definitive; Attitude Tile Header Record Contents 


Item 

Label 

Location, 
Byte Number 

Data Type 
and Length, 
Bytes 

1 

DzVIA BASE IDENTIFICATION* Index number 
referring to data base being used. 
Increments by one each time data base 
modif i ed 

1 

1*4 

2 

OUTPUT RKPETITION PERIOD. Seconds 
between attitude points 

5 

1*4 

3 

MAX PITCH. Maximum value of pitch 
contained in file in deg 

9 

R*4 

4 

MEAN PITCH. Mean value of pitch data 
points contained in file in deg 

13 

R*4 

5 

MIN PITCH. Minimum value of pitch 
contained in file in deg 

17 

R*4 

6 

MAX ROLL (deg) 

21 

R*4 

7 

MEAN ROLL (deg) 

25 

R*4 

8 

MIN ROLL (deg) 

29 

R*4 

9 

MAX YAW (deg) 

33 

R*4 

10 

MEAN YAW (deg) 

37 

R*4 

11 

MIN YAW (deg) 

41 

R*4 

12 

PERCENT PITCH DATA FLAGGED FOR 
QUESTIONABLE DATA QUALITY 

45 

R*4 

13 

PERCENT ROLL DATA FLAGGED FOR 
QUESTIONABLE DATA QUALITY 

49 

R*4 

14 

PERCENT YAW DATA FLAGGED FOR 
QUESTIONABLE DATA QUALITY 

53 

R*4 

15 

Start time of FILE (YYDDDHffllM) 

57 

1*4 

16 

End time of FILE (YYDDDHHMM) 

61 

1*4 
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Dofinitive Attitude File 
(('i.nt inuat i on 1) 


RcH'ord Contents 


Tahlo 2-4. 


Dofiuitivo Atlitudo Filo Uoador Record Conteut:s 
(Continuation 2) 


Data Type 
Location, and Length, 

Item Label Byte Number Bytea 


3b Constant a priori predicted yaw 

549 

R*4 

37-40 Constant amplitudes of Fourier series for 
a fuuorl yaw prediction Yp(Kp 1=1,2, 3, 4) 

553 

R*4 

41-44 Constant plnises of Fourier series for 

prediction Yp( i=l ,2,3,4) 

569 

R*4 

4b Fundamental frequency of Fourier series for 

a priori yaw prediction = 2n/(orbital 

period) ) 

585 

R*4 

4b Scaling factor of cross-correlation of 

589 

R*4 


errors in and in interpolation 

estimate (P" 

I IN 


47 Time constant of cross-correlation of 

errors in Y.^ and (t^^) 

593 

R*4 

48-50 Magnetometer bias field calibration 
c oe f f 1 c 1 ent s (a , b , c ) 

597 

R*4 

51-59 Nine elements of magnetometer calibration 
matrix, [M] 

609 

R*4 

..;0-71 Sun sensor alignment Euler angles 
(0.,O. ,il-.,l-l,2,3,4) 

645 

R*4 

72-143 Calibration parameters for sun sensors 
(a|,B|,J= 1,2,...9; i=l,2,3,4) 

693 

R*4 

144-306 Pitch and roll biases and telemetry 

98] - 1632 

R*4 


conversion factors plus telemetry 
conversion factors for magnetometer 
da ta 
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Table 2-4. 


Definitive Attitude File Header Record Contents 
(Continuation 3) 


1 1 em 

Label 


Location, 
Byte Number 

Data Type 
and Length, 
Bytes 

307 

Proportional contribution to 

for yaw 

1633 

R*4 


cstimace from roll (p ) 

■j 



308-427 

Various factors used in data 

processing 

1637-3476 

/ R*4 





( 1*4 
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Table 2-5. Delinitive ALLitude File Data Record Contents 


Location, 
Byte No 

Label 

I tern 
Length, 
Bytes 

Data Type 

1-24 

Logical record 1 

24 


1-2 

Day of year 

2 

1*2 

3“b 

Millisecond of day 

4 

1*4 

7 

Pitch data quality flag (See Table 2-6 
for a detailed description) 

1 

1*1 

8 

Roll data quality flag - same as for 
pitch 

1 

1*1 

9 

Yaw data quality flag - same as for 
p itch 

1 

1*1 

10 

Yaw/pitch error correlation coefficient^ 

1 

1*1 

11 

Yaw/roll error correlation coefficient^ 

1 

1*1 

12 

Yaw determination accuracy (deg, la)^ 

1 

1*1 

13-16 

Pitch, deg 

4 

R*4 

17-20 

Roll, deg 

4 

R*4 

21-24 

Y^lw, deg 

4 

R*4 

25-48 

Logical record 2 - same as bytes 1-24 
for second point on file 

24 


3433-3456 

Logical record 144 

24 


3457-3476 

Filler (all zeros) 

20 


^Data Number: 0 = ~1,00, 255 = +1*00 
^Data Number: 0 = 0.000, 255 = 0.500 
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TabU' 2-(). Description of the Attitude Data Frame 
Quality Indicator Byte 


Bit 

Number 

DesoripLioa 

Value 
of Bit(s) 

1-B 

Number of data points divided by 2 that 
were used to obtain pre-averaged value 

0-30 

() 

Attitude angle was measured 

0 


Attitude angle was interpolated 

1 

7 

Data was within "No" criteria during 
smoothing process 

0 


Data was outside criteria during 

smoothing process 

1 

8 

Nominal mode was used to determine 
attitude angle 

0 


Non--nominal mode was used to determine 
attitude angle 

1 
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SECTION III 


CHARACVI’KRI/^TION OF SKA SAT ATTITUDE HISTORY 


A. SUFmRY 

Tlio Soasat mission can be gi'oupod into tho five time periods described in 
the io 1 lowing paragraphs. 


1. Initial Acquisition of Orbital Mode (9 Days) 

The period from day 178 (Launch) through day 187 was devoted to acquiring 
a nadix'-referenced orbital attitude under control of the OACS, and to Initial 
checkout of the scientific sensors. The events of this period are described in 
Reference 3-1 aixd Section VI of Volume II of this report. By the beginning of 
day 188, all sensors had been turned oix and were operating together, and the 
spacecraft was stably controlled in attitude using only the right Scanwheel 
signal processor. Because attitude data from this period was not part of the 
scientific data record, no corrections xjere applied in ground processing. 


2. First Quiescent Data Acquisition Period (32-1/2 Days) 

From day 188 through the first half of day 220, no horizon sensor inter- 
ference occurred. The calibration of science sensors and acquisition of data 
were pursued vigorously. OACS parameter trimming was completed on day 194, 
which reduced yaw excursions by a factor of 4 (approximately). All attitude 
data from day 188 on was fully corrected in ground processing. 


3. OACS Mode Research (5 Days) 

On day 220, the sun interference with horizon seixsors resumed, as described 
in Reference 3-2. Until day 225, this period experienced a variety of responses 
to the horizon sensor difficulty, culminating in the use of Mode 5. Spacecraft 
attitude was coxxsiderably disturbed by all attempts to use the left signal pro- 
cessor. Mode 5 involved disconnection of the roll signal from the control loop, 
which proved the least disturbing method of coping with the interference. 


4. Cruise in Mode 5 (30 Days) 

From day 225 through day 255, Mode 5 was used to avoid horizon sensor 
interference. Pitch and roll were generally well-behaved, while yaw was per- 
turbed every revolution by Mode 5 to excursions of 1-3 deg, mostly over the 
southei-n hemisphere of the Earth. All orbit adjustments were achieved during 
this period, which created different patterns of attitude perturbations on the 
5 days affected by these maneuvers. Yaw attitude data gaps began to appear on 
day 226 and continued (as shown in Figure 2-12) to the end of the. mission. 
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These data gaps were interpolated with the. algorithm previously discussed. Sun 
sensor head 2 began to provide sun position data at about day 236 and by day 243 
was the sole source lU' such data (head 1 gave data until day 243). 


5. Second Quiescent Data Acquisition Period (27 Days) 

On day 255, when all chance of horizon sensor interference had apparently 
passed, the use of Mode 5 was discontinued. Normal OACS operation, on the right 
signal processor only, then resumed. Pitch and roll behavior was similar to 
that in the first quiescent period after OACS trimming (days 195-220). Yaw 
behavior was also apparently similar, but directly visible only 900-1000 s each 
revolution from sun sensor head 2 until the last week of the active mission, 
when head 3 began providing a similar span of data in another phase of each 
revolution. Most of the yaw behavior was, therefore, inferred from the partial 
data by use of the yaw interpolation algorithm. 


B. INn.IGHT ALKJNMENT CALIBRATION OF SUN SENSOR HEAD 2 
1 , Summary 

In the last week of August 1978, the sun traveled the region in which the 
FOVs of sun sensor heads 1 and 2 overlapped. Figure 3-1 s lows a sample of this 
data, with head 2 data mapped into the boresight referem frame of head 1 for 
direct comparison. Head 2 reads 0.14 deg less cross-com and 0.11 deg less cone 
than heaid 1 for the same sun position. This was significantly larger than could 
be expected from the nominal error budgets, and precipitated a detailed investi- 
gation. No conclusion could be dra\m from the overlap data alone, since it was 
obtained from such a limited area (one corner) of each head's FOV; the data 
allowed for only a two degree-of-freedor.i (DOF) solution to what was really a 
three DOF problem. 

The sun angle difference between heads remained relatively constant during 
the approximately 1-week overlap period; therefore, its cause was likely to be 
some combination of fixed changes and errors in calibration and alignment, of 
either or both heads. The largest calibration residual measured at the manufac- 
turer (Adcole) for any Seasat sun sensor was about 0.025 deg in either axis 
(cone or cross-cone) ; a worst-case combination of such errors could yield about 
0.05-deg relative difference, clearly not the answer. The calibration history 
over the I'ight acceptance test sequence shows shifts of a few hundredths of a 
degree, on the same order as the calibration residual Itself. Tolerance buildup 
could not furnish a logical explanation for the observed differences. The case 
for a fixed change or error in the alignment of head 2 was investigated and found 
compatible with observed phenomena. Though no cause for such misalignment could 
be proved, misalignment was demonstrated with confidence sufficient to warrant 
correcting data from head 2 for the calculated offset values. 
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Figure 3-1, Sun Sensor Data Dy Overlap of FOV for Heads No. 1 
and No. 2 (Sample) 
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2. Possible Causes of Data Anomaly 

There were at least two possible causes of functional misalignment for any 
sun sensor head. One was the unseating of both fine reticles, perhaps caused by 
mishandling on the ground or a flight incident. There is a precedent for this 
explanation: one fine reticle of head 4 exhibited a 0.08-deg shift over vibra- 

tion, ascribed to improper seating and corrected through rework and retest, and 
one fine reticle each in heads 1 and 3 showed smaller shifts over temperature, 
which were also improved by rework and retest. All of this occurred at the manu- 
facturer's facility and none of it involved head 2. There were no checks of cal- 
ibration performed after the units left Adcole, and, therefore, alignment stabil- 
ity was dependent on careful handling. Although the handling requirement ("... 
packaging . . . shall be such that the unit will not experience environmental con- 
ditions more severe than those specified for flight.”) would have protected 
the sun sensors if adhered to, no procedures were provided by LMSC to verify 
sensor integrity, other than the LMSC Quality Assurance (QA) stamp of approval. 
This is standard practice for low-cost projects, and leaves open the question of 
reticle shift. 

Another explanation for sun sensor misalignment is an error or change in 
the orientation of the body of the sensor head. Heads 1 and 4 were mounted 
square to the side of the Sensor Support Module (SSM) in clock and elevated 
6.5 deg above the horizontal in cone, requiring a very simple bracket mounting 
structure. Heads 2 and 3 were mounted at an angle of 23.5 deg in clock with 
respect to the sides of the SSM, and at 8.0 deg and 6.5 deg, respectively, above 
horizontal in cone. These mountings required a special small platform on which 
to mount the bracket, and in the case of head 2 several shims were required to 
reach the correct orientation. The complexity of mounting for heads 2 and 3 
made them more liable to alignment problems, including inflight shifts. There 
is no evidence other than the observed data anomalies to either indicate or vin- 
dicate the mounting as the problem source, but it remains a possibility. 

Other possible causes could be offered to explain the apparent sun sensor 
misalignment. However, no project resources could be made available to check 
them out, and the fact remains that the phenomena can be accounted for function- 
ally as alignment offsets. The following sections describe the method used to 
estimate the magnitudes of these offsets and the effects on attitude data 
processing. 


3. Alignment Calibration Method 

a. Introduction . Misalignment of sun sensors can occur in three spatial 
DOFs, and may or may not be constant in time. The sun sensor FOV overlap data 
indicated a constant misalignment. However, the data was concentrated into a 
small range of sun angles, and therefore gave an accurate measure only 
two degrees of freedom, with very pocr definition of rotation about the sun line 
(the third DOF). Another kind of data was required to pin down the other DOF, 
and it was found in the comparison of pitch reconstructed from this same sun 
sensor head with pitch computed from Scanwheel data. Pitch computed from the sun 
sensor is dependent on information in the XZ plane only, to first order, and can 
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therefore be used to unambiguously relate sun sensor cone and twist errors. By 
selecting two sun sensor data sets produced by sufficiently different sun tracks, 
this pitch comparison affords a solution for the actual sun sensor cone and twis' 
angles. Sun sensor clock for head 2 relative to head 1 can then be obtained by 
reference to data from the FOV overlap of the two heads. 

The method outlined above depends on several factors. To ascribe misalign- 
ment values to head 2, the alignment of head ’ must be in flight insignificantly 
different from ground alignment. To validly use Scanwheel pitch as a basis of 
comparison, the right Scanwheel must have had a constant set of alignment errors, 
the mapping of which into pitch must have been well measured by sun sensor head 1 
through the GSFC bias determination process. 


b. Mathematical Basis 


2 



COORDINATES 


BY DEFINITION: 

CROSS-CONE =^5= tan"' 

CONE s = tan"' (Sj/Sg) 



The cone and cross-cone angles defined above are direct outputs of the sun 
sensor. Using Equation (1), a unit sun vector (§) can be constructed from sun 
sensor data: 
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To express this vector in SARA coordinates, it is necessary to use the 
Euler angles of the alignment transformation: 


Sj^^(SARA) 


cos -sin 0 

sin cos 0 


cos 0 sin 

0 1 0 


-sin 0 cos Bg^ 


cos Ygj^ -sin Yg^ 0 
sin Yg^ cos Ygj^ 0 


S(ssi) 


( 3 ) 


where ug^, Bgi* "r'Bi clock, cone, twist, respectively, and i indicates the 
sensor head from which the data was obtained. 


The reconstruction of pitch from sun sensor data uses the observation 
model (Equation 4-166 of Reference 2-2): 


-sin p' cos r' + sin r' cos p' (4) 


where. 


P' 


= true pitch = p - p ' 

D 


r’ 


- true roll = r “ r ' 

B 


p,r = pitch, roll data from Scanwheels corrected for Earth oblat- 

ness, S/C altitude, and systematic horizon radiance 
variation 


Pb*"b 


= true bias of pitch, roll 




R1 


R2 


R3 


sun position observed by sun sensor, in SARA coordinates 
(Equation 3, above) 


^0 = 


'01 


02 


'03 


sun position determined from orbit parameters and time, 
in orbital coordinates 
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Pitch is r ccons true tec, 1 as the solution (p) of Kciuation (4) for . ' , r' , 

Sr, and Sq ^iven time. To first order the sensitivity of this determination 

of pitch to the available sources of error is as follows (see Reference 3-3): 


C.) = B “ -cos( <jj)/cos 

ijj/cos 


"4 
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-cos sin( i^^- .jj) 


ss 
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2 2 , 2 
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S B I5 
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(5) 


where S 2 and S 3 are as defined in Equation (2), and Pg) are the clock and 

cone angles of S^: 


clock = tan ^ 

“1 

cone - 3 = cos (S «) 

S Kj 


( 6 ) 


A determination of the sun sensor misalignment apparent about the pitch 
axis can be made by solution of this set of simultaneous equations: 



where Pj^j. is pitch from sun sensor data, p and pg are as defined for Equation (4), 
C‘2 and C 3 are as defined in Equation (5), the subscripts i and j refer to differ- 
ent instants of time, and Ayg) are the alignment changes in cone and twist. 

Given the relatively large stochastic errors to which p^j. and p are. subject, care 
must be taken to reduce their effects before solving Equation (7). Accurate 
solution of this equation also depends on sufficient differences between the rows 
of the C matrix, which is obtained by judicious selection of the mission times 
for which data is analysed. 
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Equation (7) solves for cone and twist independently of the sun sensor 
overlap data„ To complete the analysis and solve for clock, map the sun vector 
observed simultaneously by both sensor heads to SARA as a common reference (Equi 
tions 2 and 3) : 


s„, . [t““1 s 
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ssl J ‘ ssl 


R2 




ss2 


where T is the transformation matrix from the indicated sensor head to SARA. 
Since the two data vectors mapped to SARA represent the same sun vector: 
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+ AS 
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R2 


( 8 ) 


If there were no anomalies, ASj^2 would be zero. By considering the trans- 
formation to SARA as composed of nominal plus small angular error components, 
equations can be obtained in terms of the desired solve-for quantities: 
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(9) 


The nominal components of S]^2 (and of Sri) are obtained by the nominal 
transformation from sun sensor to SARA coordinates. Thus; 
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( 10 ) 
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Taking the difference betwetm the nominal and total (Equation 9 ) transformations 
for sun sensor head 2 gives *he difference matrix, which is, to first order 
(using yb “ 90 ”) : 


„SARA 
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( 11 ' 


Rearranging Equations ( 10 ) and ( 11 ): 
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where Si, S2- S3 are components of the unit sun vector observed by sun sensor 
head 2 . Eliminating the dependence of the first and second rows: 
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This is easily solved as: 
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A S3 + AY3 (S^ sin 3^2) 
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(14 



Therefore, cone and twist error in sensor 2 alignment can be solved from 
the sun sensor FOV overlap data as a function of clock angle (a|^2)* Solution of 
cone and twist from pitch data is independent of clock. If sensor 2 was mis- 
aligned, the two solutions will intersect at the true clock angle of the sensor 
head 2 boresight. 


c. Analysis of Data . Figures 3-2 and 3-3 show samples of pitch data 
from two revolutions on each of two vjidely separated days. The days were chosen 
to give maximum differences in sensitivity of reconstructed pitch to sun sensor 
alignment errors. Pitch is sVKn>m in Figures 3-2 and 3-3 both as reconstructed 
from sun sensor data and as determined from Scanwheel data. The latter is fully 
corrected except for the +0.11-deg bias measured earlier with sun sensor head 1; 
correction of Scanwheel data for this value of bias would widen the apparent 
separation of the two determinations of pitch. None of the data is filtered, 
but a smooth line was drawn by the "educated eyeball" method to facilitate 
comparison. 

The scatter of points for reconstructed pitch in Figure 3-2 is due to roll 
TR noise and calibration residual error in the sun sensor cone axis; in Fig- 
ure 3-2, Scanwheel pitch is scattered primarily by pitch IR noise. Between 
day 238 and day 282, the sun beta angle changed from about +51 deg to about 
-2 deg. Since the sensitivity of reconstructed pitch to roll determination 
error is approximately the tangent of the sun clock angle, the sensitivity on 
day 238 was about 1.2 to 1.8, while on day 282 the sensitivity was about 0,03 
to 0*04. Therefore, in Figure 3-2, the reconstructed pitch shows larger effects 
of IR noise than Scanwheel pitch, where in Figure 3-3 the reconstructed pitch 
shows only sun sensor calibration residual error and no IR noise. 

Figure 3-1 shows a sample of the sun sensor FOV overlap data. The compari- 
son of cone and cross-cone angles between sensor heads is subject only to mis- 
alignment and residual calibration errors. To facilitate comparison between 
the two measurements of the same sun track, the cone and cross-cone data from 
sensor 2 were mapped through nominal alignment angles to the reference frame 
of sensor 1. From this data alone, the 0.14-deg difference between cross-cone 
and the 0.11-deg difference in cone could be ascribed to raisalignments of those 
amounts in the fine reticles of one (or both) sensors, following the hypothesis 
of reticle shift. However, if reticle errors are mapped to their effect on 
reconstructed pitch, they fall far short (by a factor of 2-3) of explaining the. 
differences between the two pitch determinations. For this reason, and because 
the alternate (misalignment) hypothesis fits all the data quite well, the reticle 
shift hypothesis can be discarded. 

Figure 3-4 presents individual solutions to Equations (7) and (14), based 
on the premise of constant sensor 2 misalignment and a variety of data fx'om Fig- 
ures 3-1, 3-2, and 3-3. Each line represents the locus of cone and twist solu- 
tions available from either a single pitch comparison or a comparison of cone 
and cross-cone pairs for a single point in time. For each line, the actual 
position of the sun at that time was taken into account. Each line of overlap 
data solutions represents a range of sensor 2 clock alignment values, as indi- 
cated; the spread between lines (±0.01 deg) is due to sun sense lesJdual cali- 
bration error. The spread in the pitch comparison lines (±0«07 Ueg) is due 
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from Scanwhool data. Tho orlp.in Is marked "pre-launch" to shoX'? that 
soxipht are in terms ixf deviations from nominal alip.nmeixt. 


for day 238) 
the solutions 


rip.ure 3-5 shows the least-scpiared-error solution lines obtained from the 
data of I'lpurc 3-4. The sinp,le line for each day (238 and 282) represents the 
combination of both revolutions of data on that day. The intersection of those 
two lines is the solution to Equation (7) for cone and twist error. Not slunm in 
Figure 3-5 is a comparison of reconst I'uc ted pitch to Soanwheel pitch on days 261 
and 265; near these days, the sensitivity to twist was hw, but the data tended 
to confirm cone error of approximately -0.15 deg. (See Section 3. 1.2. 2 of Refer- 
ence 2-9.) Note that the solution locus from day 241 overlap data agrees within 
0.01 deg with the solution from pitch data at a clock alignment of 23.265 deg; 
use of the nex^ sensor 2 alignment values in GSFC processing of day 238 overlap 
data (same FOV area as day 241) shoxjed that the data anomaly had been rectified. 
Therefore, data from three different days and three different areas of sun sensor 
2 FOV all were explained satisfactorily by a constant alignment that differed 
sidnst antial ly from the groxmd-measured values. This is the simplest explanation 
consistent with all of the data, and although no mechanism or scenario can be 
offered to account for its genesis, the misalignment of sun sensor head 2 is the 
most logical cause of the observed data anomalies. The final estimates from this 
Inflight calibration are compared with groimd-cal ib rated values in Table 3-1. 


4. Effects on Definitive Attitude File 

The disagreement between sensors 1 and 2 was discovered in Septenher 1978. 
The spacecraft ceased operations in early October. Immediately after that event, 
GSFC processing of the d.ata x^as halted on project request to orgnixiue the 
processixxg of extaxxt flight data tox«ird maxlmxim sciexxtifie retxirn. Detailed 
axxalysis of the sxxix sexisor probltxm x<ras delayed xintil GSFC processing coxxld be 
resximed. Wxexx processixxg resxxmed in early November, emphasis Xixas placed oxx the. 
period betx^exnx day 256 (13 September 1978) and the exxd of the mission (day 283, 

10 October 1978). All DAFs for that period \<ro.vo completed before the results of 
the sxxn sensor misal igximent axxalysis coxxld be implemented. All yxw measxirements 
durixxg this period x^ere obtaixxed from sxixx sexxsor 2; therefore, all DAFs for this 
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Sun Sensor No. 2 Misalignment 
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period have an additional yaw error source due to that sun sensor misal ij^;ninent . 
The project decided not to request reprocessing of the 27 files already produced, 
but to change the data base for processing of earlier data. (The data from day 
220 to day 255 was next on the DAF processing schedule.) Therefore, becuuise of 
the reorganized DAF production sequence, data from sun sensor 2 early in the 
mission was processed with ground-calibrated alignment values. 

Since the differences between actual and nominal alignments have been 
determined for sun sensor 2, and the methods ChSFC used for determining yaw are 
known, the effects of using incorrect alignments are deterministic. Figure 3-6 
shows these effects for four representative days in the affected time period in 
the sense of the correction factor to be added to the DAF yaw value. The sun 
vector was within the sun sensor FOV from approximately 4100 s to 5100 s after 
ascending nod?, in each revolution (actual times are slun^7n in Figure 2-12). Dur- 
ing this time, the effect of misalignment is the mapping of the alignment error 
into yaw as a function of sun position in spacecraft coordinates (described in 
Appendix B of Reference 3-3). For those times during a revolution when the sun 
was out of the sun sensor FOV, yaw was interpolated with the algorithm described 
in Subsection lT-B-3. With this method, the effect of an error in yaw observa- 
tion falls off exponentially with time, as shown in Figure 3-6. 

The drop tt? zero correction in the north polar regicm for days 274 and 282 
does not represent an increase in accuracy. Starting on day 272 (see Figure 2-12 
for exact time), sun sensor 3 could sometimes see the sun. During those parts 
of each revolution, yaw was determined from sun sensor head 3, and, therefore, 
was not affected by sun sensor head 2 misalignment. However, analysis of pitch 
reconstructed from sun sensor head 3 data during this period showed approximately 
0.1 deg difference from Scanwheel pitch (see Refei*ence 2-9, Section 3. 1.2. 3). 

From previous arguments presented for head 2 anomalies, it would appear that head 
3 was also misaligned, at least in either or both of cone and tx^rist. Since there 
was no corresponding sun vector data from sun sensor FOV overlap during this 
period, and reconstructed pitch sensitivity to cone and tx^;ist errors remained 
relatively constant during the 9 days that head 3 was producing data, there was 
no x^^ay to separate cone from twist error, and no information at all about 
possible clock error. Reconstructed roll would have been useful in this analysis 
to separate cone from twist, but since the sun was near the orbital plane, roll 
from sun sensor data was extremely sensitive to all sun sensor errors and could 
not be trusted. In summary, no yax^^ correction could be calculated for sun sensor 
3, and the anomaly observed in data from this head could be handled only by 
increasing the allocated uncertainty for yaw determined from that source. Atti- 
tude determination uncertainties will be discussed later in this report. 
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ERECTION 





c. 


ATTITUDE HISTORY 


1. Method of Attitude History Characterization 

All of the attitude history from days 188 through 283 (EOM) is presented 
rev-by-rev in an unpublished supplement to this report.* To facilitate compari- 
sons within and between eras of the mission, a method was devised to characterize 
attitude history over spans of time, relatively longer than one revolution. This 
method was also instrumental in evaluation of yaw interpolation accuracy (to be 
discussed later) . 

Figures 3-7 and 3-8 display (superimposed) all of the attitude history 
(14 revs) of day 206. The general trend of the data can be seen as a function 
of time from ascending node, as can the deviations from trend. Figures 3-9 and 
3-10 display the mean (p) , the mean plus and minus the standard deviation (u ± o) , 
and the outer bounds of the attitude behavior of each axis. The bounds are the 
outer envelopes from Figures 3-7 and 3-8; the mean and standard deviation was 
computed at intervals of 25 s from time of ascending node. It should be noted 
here that this a is nut AD uncertainty; rather, it is a measure of the deviation 
of the actual attitude from the daily mean attitude at the given time from ascend- 
ing node. The accuracy to which these quantities are known is a separate issue 
to be addressed later in this report. 

It can be seen from these illustrations that p and a describe a kind of 
inner bound of attitude behavior, and that the. largest deviations from the mean 
are on the order of 2o and sometimes 3o. Figure 3-11 (from Reference 3-4) shows 
the geocentric representation of zero control error; since pitch and roll were 
controlled to a geodetic reference, zero control error in that frame maps to the 
non-zero function shown (perfect yaw would be always zero in both frames). The 
degree to which pitch and roll attitude behavior on day 206 met the ideal can be 
seen by comparing the characteristics showi in Figures 3-9 and 3-11. 

Figures 3-12 through 3-15 illustrate day 249 in the sa . -ay as day 206 was 

shovm in previous illustrations. (Note that the scale is t . < i for Figures 3-7 

through 3-13, while the scale of Figures 3-14 and 3-15 is r a : ; ; 'cause of the 

larger maximum excursions.) Although the method of represei.t.; '■ r. .s the same, 
the characteristic behavior is different. The biggest differu. due to the 

attitude excursions induced by Mode 5 near 3800 s, which dominates oehavior for 
the following 2000 s or so. Another difference is the brief lun izon sensor 
interference effect near 1900 s, caused by the sun passing through the Scanwheel 
active FOV. Mode 5 was used to protect against this kind of interference from 
day 248 on. Yaw was directly observable for only the 1150-s period indicated, 
and was interpolated (as previously described) for the rest of each revolution. 

The relatively narrow range of deviations in the interpolated span is due to the 
limitations of the interpolation process: the smaller the deviations from the 
mean, the less information was used in forming the estimate. 


*In cust-ody of Dr. Hiroshi Ohtakay, Guidance and Control Section (343), Jet 
Propulsion Laboratory, 4800 Oak Cnive Drive, Pasadena, CA 91101. 
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Figure 3—10* "faw Envelopes for Day 206 






Figure 3-12. Composite Pitch and Roll for Day 249 




Figure 3-13. Pitch and Roll Envelopes for Day 2.49 
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Figure 3-15. Yaw Envelopes for Day 249 


Pi^!;uros 3~16 and aro lust oj^.rains i\>r days 206 and 240 of at t Undo 

deviation from the Bioan taken every 200 s for evei*y revolution represented by 
the previous 11 histrations. Kaoh histoyj^am represents a little over 400 points, 
taken from every phase of the attitude functions. All the histop.rams fm* day 
240 show higher central peaks and longer tails than the day 206 histograms, the 
result of Mode 5 disturhances. The yaw function on day 240 is further distorted 
by the effect of yaw interpolat ion ovor nuist of each revolution; however, since 
tl\o maximum excursion region was directly observed and this region tended to 
have the highest deviations from the mt'an, the tails of the frequency distribu- 
tion are relatively undistorted. Comparing the yaw tails of days 240 and 206 
shows that Mode 5 Induced greater yaw variability. Pitch and roll under Mode 5 
show not only longer tails, but also lilgher central peaks, which reflect the 
regularity of response to the Modi' 6 coimnands. 

Figure ‘3-lB shows the probability with which an attitude deviation in a 
given axis at any time would exceed a given multiple of the overall standard 
deviation for that axis. The smooth curves sho\^m for each axis on days 206 and 
24^) were estimated graphically from the previous hlstngn'Uii data. r The same 
fiinctlon for the Gaussian vlist r ibution is included for comparison.^'^Tl^'^^Iccuracy 
of the data functions in the high-sigma, low-probability region is limited by 
the small number of samples in this region. However, there are sufficient data 
to permit the following conclusions: 

(1) For any axis at any lime, there is approximately a 7 percent chance 
of an attitude excursion from the mean function of more than 
±1.63c (for a normal distribution there wmild be 10 percent chaiu’c) . 

(2) For any axis at any time, the fiftieth ]U'rcentile of attitude 
deviations fr»nn the mean function is at abmit 0.40v' (0.67c for a 
(;aussian) . 

O) There are several times more chances for a large deviation from the 
mean function under Mode 3 conditions than under usual conditions. 

(4) For any axis at any time, the ± e envelope contains abmit 80 
percent of attitude behavior. 

Other conclusions could be dra\ai from Figure 3-18 with less confidence, but 
the above 4 conclusions are sufficient to put the \\ ± o curves into perspective. 
To the degree described here, these curves can raid will be termed "cliarac terlstic 
attitude functions." As such, they will be used to characterize spans of missloTi 
time within which conditions remained essentially constant. 


2. First Quiescent Data Acipusition Period (Days 188-220) 

This period is logically sepai'able into two subperiods, which are divided 
by the completion of OACS trimming on day 104. 
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a. OACS Trinuning Period (Days 188-194) . Figure 3-19 depicts the average 
attitude functions of time for each axis on each day of this period. Also shown 
are the contours of extreme ]i ± a for this span of time. Note that pitch remained 
the same but roll and yaw changed significantly between days 191 and 192, when 
the first changes were made in OACS parameters. As an additional aid in evalua- 
tion of these functions. Figure 3-20 shows the pitch and roll average functions 
for day 188 as uncorrected, fully corrected geodetic, and fully corrected geo- 
centric (same as in Figure 3-19) . The uncorrected function is the spacecraft 
control system view of its own behavior, in a reference system which is geodetic 
but with systematic errors. The corrected geodetic function is the actual average 
performance relative to local vertical for that day. 


b. Trimmed Quiescent Cruise Period (Days 195-220) . Figure 3-21 shows 
the characteristic attitude functions for the quiescent period after the last 
OACS parameter change of day 194 and before the resumption of Scanwheel inter- 
ference on day 220. Pitch in this period is indistinguishable from that of the 
earlier period, and remains essentially the same throughout this 25 days of cruise. 
This result was to be expected, since the previous parameter trimming was directed 
toward roll and yaw control and the only variable source of systematic error 
(horizon radiance variation) changed little on such a time scale. The roll 
characteristic fvmction can be seen to change between the beginning and end 
of this period, due to the relatively rapid shift of apogee and perigee lati- 
tudes over time and the resultant phase shift of the roll altitude correction. 

Since the Seasat orbit nearly repeated itself in longitude every 3 days, the. 
characteristic functions were plotted in groups of three consecutive cays. The 
effects of longitude variation can be seen within these triplets for roll and 
yaw. Roll/yaw control coupling is evident in the variation of the yaw charac- 
teristic, which follows the variation of the roll characteristic. The regularity 
of yaw over this period plus the roll /yaw control coupling are the basis for yaw 
interpolation in the period after day 255. (This will be discussed later.) 


3. OACS Mode Research Period (Days 220-225) 

At -approximately 08:52:00 on day 220, at the end of rev 605, sun inter- 
ference with the right Scanwheel resumed. This precipitated a period of experi- 
mentation with OACS modes of operation to find the most acceptable means of 
minimizing the effects of this interference. A method frequently used was a 
switch to the left signal processor during the time that interference was expected 
(see Figure 2-11) for the right Scanwheel, in the often frustrated hope that 
the left Scanwheel would not receive interference at the same time as the right 
Scanwheel. This mode of interference avoidance was designated Mode 12. Also 
attempted xo^ere all combinations of primary and secondary Control Logic Assembly 
(CLA) power supplies. During revs 626 to 628 on day 221, both right and left 
Scanx^heels were used simultaneously in the intended normal mode of operation; 
Figure 3-22 shows the effect of interference in this mode. 
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-20. Corrected Versus Uncorrected Pitch and Roll (Mean for Day 188) 




Figure 3 - 21 . One-Sigma Attitude Envelopes for Days 195-219 



Figure 3-21. One-Sigma Attitude Envelopes for Days 195-219 (Continuation 1) 





Figure 3-22. Examples of Scanwheel Interference (Days 221 and 223) 
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^ure 3-23. Examples of OACS Mode 12 (Day 220, Revs 608-611) 
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Table 3-3. OACS Mode Changes - Days 220 Through 225 (Continuation 1) 


Day 

(1978) 

Event (GMT), 
H:niln:s 

Mode Status Duration 
(Seconds From Ascending 


Mode Select 

Statufj 

• Pitch and Roll 
Bias Set 

Node) 

CLA PoT»er Signal 

Roll Attitude 


From (Rev) 

To (Rev) 

Supply 

Processor 

Signal 


221 

09:51:00 

5251 (620) 

5851 (620) 

2 

Right 

Disconnected 

2 

221 

10:01:00 

5851 (620) 

5213 (621) 

2 

Right 

Connected 

2 

221 

11:31:00 

5213 (621) 

256 (622) 

0 

Right 

Disconnected 

2 

221 

11:49:00 

256 (622) 

5296 (622) 

2 

Right 

Connected 

2 

221 

13:13:00 

5296 (622) 

578 (623) 

2 

Right 

Disconnected 

2 

221 

13:35:00 

578 (623) 

165 (626) 

2 

Right 

Connected 

2 

221 

18:30:01 

166 (626) 

3690 (628) 

2 

Both 

Connected 

5 

221 

22:50:00 

3690 (628) 

1852 (629) 

2 

Right 

Connected 

2 

222-224 

00:00:00 

1852 (629) 

1432 (672) 

2 

Right 

Connected 

2 

225 

00:00:00 

1432 (672) 

4648 (675) 

2 

Right 

Connected 

2 

225 

05:55:30 

4649 (675) 

5549 (675) 

1 & 2 

Left 

Connected 

4 

225 

06:10:31 

5550 (675) 

4633 (676) 

1 & 2 

Right 

Connected 

3 

225 

07:35:53 

4634 (676) 

5534 (676) 

1 & 2 

Left 

Connected 

4 

225 

07:50:54 

5535 (676) 

4619 (677) 

1 & 2 

Right 

Connected 

3 

225 

09:16:16 

4620 (677) 

5520 (677) 

1 & 2 

Left 

Connected 

4 

225 

09:31:17 

5521 (677) 

4604 (678) 

1 & 2 

Right 

Connected 

3 

225 

10:56:39 

4605 (678) 

5505 (678) 

1 & 2 

Left 

Connected 

4 

225 

11:11:40 

5506 (678) 

1423 (681) 

1 & 2 

Right 

Connected 

3 

225 

15:05:30 

1423 (681) 

4835 (682) 

2 

Right 

Connected 

2 

225 

17:43:00 

4835 (682) 

5315 (682) 

2 

Right 

Disconnected 

2 

225 

17:51:00 

5315 (682) 

4826 (683) 

2 

Right 

Connected 

2 

225 

19:23:28 

4826 (683) 

5306 (683) 

2 

Right 

Disconnected 

2 

225 

19:31:28 

5306 (o83) 

4816 (684) 

2 

Right 

Connected " 

2 

225 

21:03:56 

4816 (684) 

5296 (684) 

2 

Right 

Disconnected 

2 

225 

21:11:56 

5296 (684) 

4806 (685) 

2 

Right 

Connected 

2 

225 

22:44:24 

4806 (685) 

5286 (685) 

2 

Right 

Disconnected 

2 

225 

22:52:24 

5286 (685) 


2 

Right 

Connected 

. 2 
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The yaw excursion was precession due to the roll torque acting on the pitch 
momentum bias. The pitch interference effect disappeared under Mode 5 (see 
Section 2.4 of Reference 3-5). The fact that the interference phenomenon 
moderated so strongly when the roll control signal was disconnected is a clue 
to the cause of the interference, but no resources have been applied to analyze 
this. Mode 5 was used to avoid the interference phenomenon from day 225 through 
day 255, when the geometry of the sun track finally precluded any intersection 
with the Scanwheel FOV. 


4. Period of Cruise in Mode 5 (Days 225-255) 

Mode 5 was used on a regular basis to moderate the effects of Scanwheel 
interference until the geometry of the sun track finally precluded any intersec- 
tion with the Scanwheel FOV on day 255. The time points of roll control discon- 
nect and reconnect were adjusted regularly to track the predicted times of such 
interference. Figure 3-25 shows the characteristic functions and the schedule 
for these mode switches over the early part of this period. The time span of 
disconnected roll control was 8 min from the latter part of day 225 until the 
end of day 226, and was 6 min thereafter. The phasing of Mode 5 was as shown 
in Figure 3-25. The experimentation to find the most beneficial phasing is 
evident both in the schedule and in the results (decreased excursions in roll 
and yaw) though the results are also aided by changes in sun track geometr-y. 

Figure 3-26 shows the characteristic functions and Mode 5 schedule for the 
latter part of this period. On day 247 the sun track began to intersect the 
active FOV of t?.. jcanwheel (near the Earth horizon), producing both pitch and 
roll disturbances near both leading and trailing edge horizon crossings. This 
phenomenon was expected from pre-launch analyses, and it lasted for 8 days. 
Beginning on day 248, Mode 5 was used for a 6-min span twice each revolution, 
as shown in the schedule of Figure 3-26. Interference in the active FOV had an 
obviously stronger effect on Scanwheel performance than interference in the 
blanked FOV, as can be seen by comparing Figures 3-25 and 3-26. Pitch by then 
had a marked response to both FOV crossings, roll response was about the same as 
previously, and yaw response was stronger. Because there was no sun sensor 
coverage (see Figure 2-12) for the northern FOV crossing, the yaw response to 
Mode 5 for that interval is not directly determinable. Figure 3-26 shows the 
results of the interpolation algorithm, which is of too low an order to accurately 
compute such a transient. The tightness of the la boundary in this region is 
further evidence of the algorithm's low order. The southern FOV crossing was 
largely observed by the sun sensor and, therefore, is accurately portrayed. The 
roll response to the northern Mode 5 application is rather weaker than to the 
southern. An analysis of wheel speeds, if it were available, would probably show 
that Scanwheel speeds were near the bias value at the. northerly roll disconnect 
and relatively far from bias at the other disconnect. 

Figure 3-27 gives two samples of attitude behavior under control of ascent 
stage gyros with torque supplied by thrusters. The examples chosen were the 
calibration burns on days 227 and 235, and the plots are phased relative to 
start of the burn. Each burn was 60-s long, and the two burns were in opposite 
directions. This fact can be seen in the attitude response to thruster 
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Figure 3-25. One-Sigma Attitude Envelopes for Days 22h-23h 
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Figure 3-26. One-Sigma Attitude Envelopes for Days 242-251 








miGalignraent with respect to spacecraft center-of-gravity . Note the -2 deg 
offset of yawj caused by gyto misalignment. The transient observable at bum 
plus 2400 s on day 227 coincides with no recorded command. The switch from 
thruster/gyro to reaction wheel (normal mode) control occurred a little after 
burn plus 3000 s on both days, and resulted in a smooth changeover. 


5. Second Quiescent Data Acquisition Period (Days 256-283) 

Figure 3-28 shows the characteristic functions for the latter part of the 
mission. It can be termed quiescent because it includes no sun interference, no 
Mode 5, and no maneuvers. The operational mode was the same as in the first such 
period (days 195-219) . Comparison of Figures 3-21 and 3-28 shows many similari- 
ties and some interesting differences. Pitch is almost identical; roll differs 
primarily in the phase of the altitude correction; yaw in the later period is 
mostly interpolated based on the average behavior in the earlier period, and 
looks similar for that reason. Yaw was directly observed (see Figure 2-12) by 
sun sensor iiead 2 near the south pole during each rev in this period, and by sun 
sensor head 3 near the north pole beginning around day 274. Note that the yaw 
interpolation algorithm, due to its lower order, tends to supply essentially the 
same estimated yaw function of time for all yaw data gaps. Tlie accuracy with 
which thiS estimate represents true yaw behavior will be examined in subsequent 
sectionL , 
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SECTION IV 


YAW AD PARAMETER ESTIMATION 


A. INTRODUCTION 

In Subsection II-B-3, the Seasat yaw interpolation method was described. 

For the algorithm to have acceptable accuracy, JPL was required to specify for 
GSFC the values and time period of applicability of the constant parameters in 
Table 2-3. The JPL analysis that estimated these parameters was based on the 
repeatability of yaw as a function of sub-spacecraft latitude and as a function 
of a simplified model of yaw/roll cross-coupling. Data from the first trimmed 
quiescent cruise period (days 195-219) were used to estimate parameters for the 
second such period (days 256-283) . Data from the early part of the Mode 5 cruise 
period were used to estimate parameters for the latter part of that period. Data 
from the first two propulsive maneuvers (days 227 and 230) were used to estimate 
parameters for all periods under RCS/gyro control. 


B. DATA PROCESSING OVERVIEW 

A set of undocumented developmental programs was used to help estimate yaw 
AD parameters from GSFC DAFs using the JPL Univac 1108 computer. The first step 
was to copy the DAF from the GSFC-supplied AOT tape to a Univac 1108-compatible 
tape; this process was performed on the JPL IBM 360 computer. Then the DAF (still 
in the IBM format) was processed by a Univac 1108 number-crunching program to 
produce (1) an edited (time and attitude values only) file in Univac 1108 "blank" 
format of all of the attitude records for the given day of data, and (2) a file, 
also in blank format, of the mean and standard deviation for each axis at 242 
equally spaced phases of the orbit (integer multiples of 25 s from ascending 
node) . The mean and standard deviation attitude functions were used to character- 
ize attitude behavior as described in Subsection Ill-C; the mean function was 
also used in another Univac 1108 program to find the best (least-squared-error) 
fit truncated Fourier series for that day of data. The parameters of several such 
(daily) Fourier fits were evaluated to determine the most appropriate fit over 
several days to estimate the parameters for the yaw interpolation algorithm. 

The edited file of all of the attitude data for a given day was processed by yet 
another Univac 1108 program to find the least-squared-error fit to the four para- 
meters (two of magnitude, two of phase)of the simplified yaw/roll coupling model. 
The parameters of several of these fits were then evaluated to determine the most 
appropriate set to use for K;^, tp^, and tj^Q in the yaw interpolation scheme. 

Because the Seasat OACS time constant (equivalent) in yaw was 3 or 4 times 
as large as the roll axis time constant, roll data had to be significantly smoothed 
before estimating the yaw/roll coupling parameters. The JPL estimation program 
used a cubic spline fit to roll on 300-s centers; GSFC AD software used a 
Chebychev polynomial instead of a spline fit but obtained similar results 
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TUo luu'oi taint ol the Fourier fit prinliet aiul ol the yaw/ri\ll eiuipliui> 
model were also etuahl lfd\ed. For the predict model, the following, variances were 
sumnu'd : 

U) Actual attitude in the refereuci' era with respi'ct to its mean 
at t i t ude 1 uuc t ion • 

Avera>\e tiueaiO luuction with respect to its Fourier approxlmat ion* 

C‘^') Overall Fourier approximat ion with respect to the averap,e function 
of available attit\ule data in the tarp.et era. 

In the latter part ol the Mode S cruise period, there were no appHcahle 
data irom which to estimate item (3), so it was estimated as an arbitrarily 
larp.e valvu,' tor interpolat ions in that period only. The variance ol the yaw/ 
roll coupling, model was estimated in a manner similar to that described for the 
predict, but itian was not included since actual, rather than averap,e, data 
was used direct Iv in estimating the parameters. The uncertainties i\st imated 
were used diri'ctly in the yaw interpolat ion algorithm. 


0. DATA PROt'KSSlNO DKTAIL 

1 . Four ier Series Approximat ion 

To obtain a predict ion function for yaw behavior, actual yaw data over a 
relerence period ol available data was characterised as a truncated thins through 
Unirth harmonic') Fourier series based on time from ascending node. Although the 
orbital period varied over the course of the mission (Figure ‘VI), this created 
no l irst'”order error for the predict function, siitce the Fourier series was 
t^xpressed in terms ol orbital c<'ntral angle. Therefore, although the attitude 
data was expressed relative to time, It was analyr,ed and character ized in a t Ime- 
noimui 1 i red st'use . 


The program YAWFIT (which obtained the least ’-squared-error lit ol the trun- 
cated seriL's to a tile ol average yaw attitude versus time) would first subtract 
the I'liect ol solar torque lor the day the data was taken and add the solar torque 
eflect lor a given target data. The solar torque model used lor thds calculation 
wa.s a simplll led version ol that published by LMSC as Reference 4-1, and the 
ellects on spacecraft attitude were I'omputed through a linearir.ed OAC'S model as 
published in Relerence 4-2. Required as inputs to the program were sun beta angle 
(Figure 4-2) and the orbital phase angle of sun nadir (Figure 4~3) for both the 
day the relerence data was taken and for the day for which the predict would be 
madt'. Tlu' solar torques on the two days were generally quite different relative 
to each other, but in practice they accounted for only a minor proportion of total 
yaw behavior. 
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a, c^u it’ina'ut IN\r iod. To pwpnw^ a pianiiot liiiu*t imi ior tho aiHunul 
quioiUHait |naiod (.ilavj^ , tlu' 'iata ol tha first .stirh pt‘riiul (days 19‘>-2I9) 

was analy:a*tl. idy^uri' ^-21 t'xampK's of tlu* avoray.o altiludt* fxmcl Ions in the 

roforonoo porhui* Ovor 80 poroont o( :.lu‘ o.vailaiilo data in this poriod was nsod 
in tho analvxsi;;. IMot.s of individual Fmirior torm oooffioionts vorsus data day 
showod soino ovolio boliavii>r ovor tho 2^)~day span, hut tho information was Insuffi- 
oloiU to lustily c'xtrapoliit ion of tlu- oyoios 1 1 ' tho larp^ta ora 30-80 days later; 
thoiadoro, tlu- nu\in value over an appaiaait half--oyole was solootod for each 
compcment tt^ form tho Fourier stories characttu' i st i c of the reference era* 

fdp^ure •4-*a c. ives two toiamjHes of the aver ay, e yaw function for a day 
in the roforonce [uu'iod vi’rsus lua h the truncated Fourier appn^xiiuatlon for that 
day and tlii* truncated Fourier seriot^ used ti> cliaractor rze the tnitire reference 
period. Tlu' offt'ct of sorites truiueititm at fourth t>r!>ital harmonic can be seen 
in the ctanparistm ta dailv avorviy,e function with dailv Fourier approximat itm; 
overall standard deviatiem t>f thi.s effect was near 0,04 dey. The differences 
between the individual Fourier appiaiximat i tuis ftu" davs within the reference period 
and the overall Fmirier series fi'r that entire period has a standard deviation of 
about 0*10 doy. If tlu‘ variatiem of actual yaw data with respc*ct tt> the average 
yavN^ fuiuH itni ft>r each day is taken itUa^ account, it is ftumd that the overall 
truncated Fourier series can idiaracteri^u^ the refertuice {^c^riod yaw behavior with 
a statistical accuraev ef about 0.2/ deg (10, 

This accuracy of cTiaracter i ;:jU ion is ttctaially a functiem in" orbital 
phase, ranging Iron: 0,11 deg to 0.40 deg (10, and is shln^m as an envelope about 
the overall Fourier function in Figure 4-3. To iiulicati' this functionOs accuracy 
in predic.ting yaw behavior for the second quiescent period (days 256-283), the 
average yaw functions fiU" several selected days within the later period are also 
plotted in Figuri' 4-5. Yaw was directly observable for only a part of each 
revolution because of sun sensor FOV limitations. The yaw plotted between 4100 
and 5100 s was corrected for the alignment errors deteniiined for sun sensor head 2 
(see Figure 3-6)- (Without this correction, these yaw^ functions would be approxi- 
mately 0.2 deg more negative; that the correction improves the comparison with 
the first quiescent period adds confidence to the aligjiment analysis.) The yaw 
plotted between 850 and 1 900 s is from sun sensor bead 3; this head was apparently 
misaligned also (see Subsection 1 I i-C) > but the range of data is insufficient to 
estimate the parameters. Therefore, the average data from this head should he 
viewed as subject to constant errors of undeterminable magnitude and direction, 
but probably within an envelope of ±(),25 deg or so. This error is, of course, in 
addition to the normal variation of yaw about its mean function. 

From the data in Figure 4-5, it can he determined that the RMS error 
between the overall Fourier F(t) and the yaw functions for days 273 and 282 is 
about 0.44 deg, if only sensor head 2 data is used for the comparison. On other 
days, this error is li^ss, but, for conservatism, the larger value was assumed. 
Taking this RMS error as one standard deviation and combining it with the overall 
op from Figure 4-5 leads to a value of 0.55 deg as the conservative lo accuracy 
of F (t) in predicting yaw for the period of days 256-282. 
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Figure 4-4. Daily Average Yaw versus Fourier Approximation - 
First Quiescent Period 
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b. Mode 5 Period . In a manner similar to that described for the 
quiescent period, a truncated Fourier series was constructed to predict y.w for 
the varied yaw data gaps that occurred over the period of days 226-255. Such 
gaps began to appear at the very beginning and rapidly became so large that only 
days 226-229 have enough data to support Fotirier analysis. Day 22V was excluded 
becauvSe several revolutions within that day were on RCS/gyro control, which has 
a quite different characteristic. The Fourier function which best fit days 226, 
228, and 229 was then adjusted to emphasize the orbital phase covering the first 
4000 s from ascending node and consequently de-omphasizlng the last 2000 s. The 
reasons for this were threefold: (1) the period between 4000 and 5100 s from 
ascending node was usually covered by sun sensor data, so interpolation was not 
usually necessary; (2) toe Mode 5 yaw excursion occurred regularly in the latter 
part of each revolution, and the truncated Fourier series cannot fit such a large 
excursion without being a poor fit to the rest of the data; and (3) the Mode 5 
excursions varied significantly in orbital phase and amplitude, which meant that 
no constant predictor function could make a good fit except during a day or two 
of this mission interval. 

Figure 4-6 shows the Fourier function chosen to represent yaw for the 
entire Mode 5 period, compared with the average yaw functions from the 3 days used 
to construct the function (Part 1 of Figure 4-6) and with the average yaw functions 
from several selected days subsequent to day 230. The average yaw functions are 
plotted only for djita from the sun sensors, so the yaw data gaps are obvious. 

Part 2 of Figure 4-6 shows clearly how different the yaw behavior was in the latter 
part of the Mode 5 period compared to the beginning. 

Figure 4-7 plots the RMS error of the selected overall Fourier func- 
tion with respect to averiige yaw as a function of time from ascending node. Note 
that the lower curve, representing the time period between days 225-235, shows 
fairly good performance in the first 4200 s of a revolutloni, and much poorer per- 
formiince in the interval containing the Mode 5 excursion. This was the intended 
result, for reasons given previously. The upper curve, representing the latter 
portion of the Mode 5 era, shows very much worse perfornuince near the Mode 5 
excursion, suggesting that performance in the orbital phase not observed by the 
sun sensor auiy also be significantly worse. 

T’r data from Figure 4-7 was used to estimate the accuracy of the 
overall Fourier function for yaw Interpolation in the Mode 5 period. Table 4-1 
lists these estimates as functions of orbital phase and mission phase. Since the 
Fourier function was used to Interpolate yaw primarily in the interval between 
0 and 4000 s from ascending node, the 0.29-deg figure was used as la for the day 
226-235 period. There was no corresponding estimate for the day 236-255 period, 
because of the lack of sun sensor data in the appropriate orbital phases, so a 
default value of 0.75 deg (la) was used. 
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Figure 4-6. Overall Fourier Function versus Daily Mean Yaw for 
Mode 5 (Continuation 1) 
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Table . Tmirier Funetion Aeeuracy for Mode 5 Period 


RMS Krror (deg) of Given Funelion per 
Orbital Phane (Seconds from 
AvScending Node) 


Fuiu’ t i im:i 1 Ri> 1 a t i onsh Ip 

Days 

226 - 23') 

Days 236 - 255 


0 - 4000 

4000 - 6000 

4000 - 5200 

Overal 1 Fourier functimi relative 
1 average vaw function 

0.16^ 

0.53 

1.27 

Actual yaw relative to average 
yaw f unct ion 

0.2 5 

0.39 

0.42 

Overal 1 Fourier function relative 
to actual yaw 

0.29 

0.68 

1.34 


‘ RG S C o n t r o 1 P e r 1 o d s ( Ma n e u ve r s ) . Yaw characteristic behavior under 

RCS/^yro control was quiescent, with small excursions relative to a mean value 
ol -2.1b dt‘p,. This large mean offset was due to an unsuspected gyro misalignment. 
Since the dvnamic behavior was apparently not a function of latitude or orbital 
phase, the Fourier function used to predict yaw for maneuver periods was simply 
the offset value, with all sinusoidal terms zeroed. The accuracy of such yaw 
estimation was calculated as 0.71 deg (Ic) by comparing the constant offset model 
with most oi the available data from RCS control periods. 


2. Roll/Yaw Coupling Mtnlel 

IkuMuse of the pitch momentum bias, roll and yaw were coupled gyroscop ical ly 
in addition to the inertial quarter-orbit coupling of a nadir-referenced space- 
craft. i’he model uhed to estimate yaw from these relationships was a linear 
combination of roll and roll rate functions: 

Y 3 - cKyr 4- tj^) -f (l“-f)KYRD R(t +tj^j)) 

Paranu^ter estimation was performed from data sets, each composed of a complete 
day's worth of data taken every 100 s. Yaw and roll data were each fit with a 
tnibic si>line, yaw on 200-s centers and roll on 300-s centers. The fit accuracy 
(hO during quiescent periods was 0.002 to 0.003 deg for yaw and about 0.02 deg 
for n^ll; for Mode 5 periods the fit accuracies were 0.01 to 0.03 deg and 0.03 
to 0.03 deg, respectively. The lower fit accuracy for roll indicated that higher 


froquoncy roll oomponontvS wore boing fiUorod out by this process, which was 
clcsixMblc since only the lower frequencies were likely to be com»laled with yaw. 
Evaluation of the spline at any time (t) produced an accurate value for position 
(Y or R) as well as rate (R) . The program used to estimate the parameters of the 
model in gain (eKyj.^ and (I-p)Kyrd) and phase (t|^ and would search for the 

local minimum in the variance between Y^, as defined above, and yaw data. For ea 
data set, three separate evaluations were made: one each for roll position (p=l) 

roll rate (e-0) , and the combination of roll position and rate. From several sue 
evaliuUions a representative set of parameters was selected, and Y*j from this 
representative set was compared again with several data sets to estimate the over 
al 1 periormance level . 


a. Quies cent Period . Data from several days in the first quiescent 
period was analyzed to produce parameters for yaw interpolation in the second 
quiescent period. Since the. orbital phase from 4100 to 5000 s after ascending 
node was normally observed by sun sensor head 2 in the latter period, yaw inter- 
polation would not have to be performed for that orbital phase. Therefore, that 
Interval was excluded from the parameter estimation process. T"he results 
suggested that tor roll position alone, the gain (Ky^) should Oe 1.48 and the 
phase {IyO about 1000 s. For roll rate alone, the gain (KyRi)) should be 7.50 
with phase (Ird) of 450 s. In combination, the value of 0.4 for the relative 
proportional factor (|. ) gave the best results. When this combination of para- 
meters was checked in the Y 3 model against actual yaw data, for days 205, 214, 
and 218 (exclusive of the 4100- to 5000-s phase), the RMS error was 0.34 deg. 
Performance was also checked against available (preliminary) data irom the second 
quiescent period. This data, from days 268, 274, and 282, had to be corrected 
for roll bias and sun sensor 2 misalignment before comparing with early data taken 
from sun sensor 1. However, the corrected data yielded RMS errors against the 
Y 3 model of 0.27 deg, which is the same result as a similar check against data 
from the same orbital phase in the first quiescent period. Therefore, it appears 
that the roll/yaw coupling model is as accurate on early quiescent period data as 
on the I a t o r da t a an d , the re. f cjr e , c an be used wit h c on f id enc e to i n t e r p o 1 a t e yaw 
dill A gaps. 


Figure 4-8 shows two revolutions of yaw data from day 218, together 
with the Y 3 function from the corresponding roll data for each revolution. The 
revolutions were chosen to illustrate relatively good and relatively poor corres- 
pondence between Y 3 and true yaw. Note that Y 3 tends to stay in phase with yciw 
but can have significant amplitude errors. Figure 4-9 makes the same kind of 
comparison in the target era: the second quiescent period. As far as can be 
determined with the limited data available, Y 3 compares to yaw equally well in 
both ( I u i e s c e n t periods. 


b. Mode 5 Perio d. Because of the rapid growth of yaw data gaps in the 
early part of this period, only days 226-229 have enough data to justify para- 
meter estimation. For the same reasons described for the quiescent periods, the 
interval between 4100 and 5000 s from ascending node was excluded from the 
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estimation process. Day 227 was also excluded since a signi^'lcant fraction of 
that day was spent on RCS/gyro control for an orbital trim mancfiver. The analysis 
showed that the best parameters for roll alone were 1.50 in gain (Kyj^) and 700 s 
in phase (t^). Ft^r roll rate alone, the best parameters were 925 in gain (KyRi)) 
and 290 s in phase (tp-)). The best combination of the two was obtained with the 
proportionality constant p=0. 2. For the day 226, 228, and 229 data, exclusive of 
the 4100- to 5000-s interval, the performance of this combination was about 
0.26 deg RMS. 

In the transmittal of these parameters from vIPL to GSFC, the Ky^D and 
values were inadvertently interchanged. All of the Mode 5 vaw data from day 
227 through day 255 was interpolated, where required, with %d= 700 and tR=925. 

The otlier parameters were as given above. Recheck of this (actual) yaw/roll model 
against the day 226, 228, and 229 data shows performance of 0.29 deg RMS rather 
than 0.26 deg RMS. Although this represents a loss of performance, it is not 
serious., and the difference is within the accuracy of the parameter estimation 
process . 


Figure 4-10 shows three assorted revolutions of yaw data from the 
early Mode 5 period along with the Y 3 functions from the corresponding roll data. 
Y 3 is shown for both the best set of parcimeters and the set actually implemented. 
The difference in performance is difficult to see except during the Mode 5 yaw 
excursion. Note that the Y 3 function has an erroneous negative transient which 
is due primarily to roll rate behavior at the Mode 5 event; the Y 3 model describes 
fairly well the steady-state roll/yaw relationship, but does a relatively poorer 
job on transients. The accuracies quoted earlier do not include this transient 
response error. 


Figure 4-11 gives two examples of yaw data compared with the r^ctual 
Y 3 function for the latter portion of the Mode 5 period. Yaw data was limited to 
the 4000-5000 s interval for which Y 3 was not optimized, and, therefore, the 
comparisons are not favorable. Also, the phase of the Mode 5 event advanced 
regularly until it could no longer be observed by the sun sensor. Table 4-2 shows 
the accuracy of Y 3 in the 4100- to 5000-s interval for both early and late portions 
of the Mode 5 period. Those accuracies reflect performance under transient condi- 
tions; most of the revolution is more nearly at steady state, and the Y 3 pariimeters 
estimated from the early data can safely be assumed to represent yaw (outside the 
transient region) as well in the later portions of the Mode 5 pericui as in the 
earlier portions . 


• RCS Control Periods (Maneuvers) . When tlie spacecraft was under 
RCS/gyro control , the misalignment of the yaw gyro caused a large fixed offset 
in yaw. Since the Y-j mcHlel had no ['jrovisions for such offsets between mean roll 
and/or roll rate and mean yaw, the Y^ model had over 2-deg RMS error during such 
periods. Thereft^re, Y 3 was not used for yaw interpol a t i on during numeuvers. 
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Table 4~2. Y.^ Accuracy (deg) at 4100-5000 s From Ascending Node 


Day of Year 

Mean Error 

Standard Deviation 

RMS Error 

226 

-0.75 

0.53 

0.92 

228 

-0.62 

0.35 

0.71 

229 

-0.67 

0.34 

0.75 

Overall 
(Early Mode 5) 

-0.68 

0.42 

0.80 

245 

-1.36 

0.69 

1.53 

254 

-0.43 

0.55 

0.70 

Overall 
(Late Mode 5) 

-0.90 

0.63 

1.09 


D. DATA BASE FOR YAW ATTITUDE DETERMINATION 

Table 4~3 lists the values of yaw interpolation parameters and sun sensor 
head 2 alignment angles in each of the five AD data bases created for Seasat. 
Also given is the range of days of data to which each data base was applied. 
Table 4~4 lists the per-revolution applicability of data bases 3, 4, and 5 on 
maneuver days. 


4-20 


r 












y-s 












in 

O 


1 




in 

rH 

o 




O 

VO 











<5* 

IM. 

(JV 




o 

o 

r»«. 

m 

o 

in 

o 


rH 

”5 

Cy 

o 

o 


1—4 

or 

O' 

o 

CM 

CM 

rH 

00 




o 



f 

o 

CM 

O' 

O 





« 

• 








in 

in 

o 

o 

f— ^ 

r*-. 

O' 

CM 

rH 

CM 

d 

d 

d 

d 

f-H 

cn 


O 

d 

d 

CO 

CO 

O' 





w 

w 





W 

1 





cc 

J 



in 



CM 

O' 

00 

1 

» 




CM 

CM 





CM 



















V 

y— » 






O 

CM 


oo 

O' 







in 

o 






m 

m 

o 




VO 

vO 

GO 

MX 

m 

m 







ret 


O' 




o 

o 

h- 

in 

o 

in 

o 

C'J 

CM 

cn 

CM 

<n 


rH 


00 

O' 

o 

CM 

CM 


00 




o 



• 

o 

CM 

ex' 






• 

* 

• 








k 

0) 

Mf 

o 

o 

O 

rH 


cn 

CM 

d 

O 

o 

O 

o 

o 

rH 

cn 

MT 

o 

d 

O 

CO 

00 

O' 

1 

w 


rH 

>w» 








1 

1 

1 

1 

00 



in 



CM 

O' 

00 


w 

u 

cu 

4J 

cu 

J-l 

td 

cu 

CO 

cd 

P3 


cd 

Q 


CO 

I 


r£J 

cd 

H 


0 

1 


<M f\J fN 


o 

o 


CO ^ 
'4- 


crt m 00 




(U 

s 

B 

o 

>-l 

U-i 

•H 

ca 


Oi S Q 

'•CNQcvJO >■>’(36 gcNn O »-< 

t2 u u 


F-< CM ro 






00 

OO 

00 





oo 

3 

3 

3 





3 

no 

no 

3 





T 








M 

i 

& 

s 

00 




3 

3 

3 

3 

3 


u 


3 

>> 


>, 

3 


Qi 








tn 



T3 

3 

3 





no 

3 

3 

3 

> 


00 


3 

4J 

4.1 

4J 

3 


(U 



U 

U 

O 

Pm 


•3 


O 

H 

•H 

•rH 





•H 

3 

T3 

3 

3 




»3 

3 

3 

3 

3 


00 


3 


U 

bl 

U 

00 

1) 

CM 

(H 

a 

CU 

a 

CJ 

0) 

•3 

00 

CL 




■H 

♦o 


3 


n 

U 

u 

3 




U 

o 

O 

0 

3 

00 



O 

Uh 


llH 

M 

(U 

■u 

CO 

4h 




a 


T3 

rt 

M 


bO bO 
di 0) 
TJ no 


s 


0) 0) -xa 


to 0) 0) 0) 


to <0 (0 CO 


(U 

(U 

bO 


cn 00 o\ 
CM O' 00 
I 


moo 


00 O 
O' ex' 
\ 


CM 

y^ 

m 

o 

CM 

y-> 

<n VO 

o 

o 

O 

CM 
y-\ 
*— 1 

o 

o 

o 

o 

■D 

o 

o 

o 

o 

3 

3 

iJ 

0 

Cli 

o 

o 

O' 

o 

in 

O 

m 

CM 

O 



o 

o 

o 











0 






o 

o — 

tn 

<r 

sr 

O 

d 

d 

d 

d 

d 

O' 

d 

d 

d 

O 

d 

o 

cn 

00 

d 

■w 


•-* 

CM 








Mf 




U 



CM 

ex' 

CTt 


t CM CM 
CD 


0} 

§ 

u 

td 

a 


td 

rH 

o 

B- 


*t3 

<i) 

>* >-• 

I § 

>> 


0 

o 

u 

tw 

§ 

f>> 


§ 

>> 


:3 

td 

g 


to to CO CO 


g 


g 


(U 

J-l 

a 


td| 


U M 
3 

o o 


0 0 0 3 

a g E ^ 

a) a) Q> CM 


-,0 O 

o 3 : ■» 

3 CM cn • 


QJ 

*3 

3 


c 

o 

CJ 


3 

o 

CJ 


3 

o 

u 


rl ri I— I OT 

0> Dt ^ & td 

I I I I £ 


CM 

tW 

O 


3 

3 

3 

t 






U 

•H 

•H 


3 

MH 


OO 


00 

3 

M 

M 

u 

3 



3 

00 

3 

3 

3 

3 

3 

O 

m 


3 

3 

3 

>V 

3 

W 

•H 


>> 

3 


3 





U 

U 


3 




iw 

bi 

ILi 

3 


c 

1 — 1 

CM 


CM 

0 

3 

3 

O 

IH 

4J 

00 


CM 



•H 

tH 

4t 

>• 


3 




3 


M 



3 

3 

3 


3 

>' 

3 

3 

MH 

3 

O 


3 

3 

3 

3 

O 

O 

0 

3 

•H 

(-> 

3 

3 

3 

3 


44 


3 

U 

3 

rC 

3 

43 

N— <• 



Pm 


3 

H 


43 



H-l 

<4H 

u 

m 

rO 

3 

u 


U 

Pm 

0 

O 

c 

pw 

•H 

W 

0 

bi 

0 

4-1 



3 


)H 


3 

O 

3 

•H 

a 

a 

3 

UH 

u 

4J 

3 

3 

3 

rH 

u 

0 

cr 

O 

3 

3 

3 

3 

3 

•H 

3 

3 

3 


O 

3 

3 

3 

3 



O 

(U 

to 

3 

X 


00 

5 


I 


3 ^ 

<S 5 


i 

£ 


4-21 


226 283 235 255 Table 4-4 


Table 4-4. Data Base Applicability for Maneuvers 


DOY 

Rev 

Event 

From 

Times 

To 

YAWINT 
Data Base 

Event 

Description 

227 

700-702 

00:00:00 

4:18:20 

3 

Normal (Mode 5) 


702 

2:30:19 



Start RCS 


703 

4:19:50 

4:29:10 

5 

Yaw Data Gap 


704 

6:00:15 

6:09:55 

5 

Yaw Data Gap 


705 

7:41:08 

7:42:08 


Burn 


705 

8:32:00 



Stop RCS 


705 

7:40:30 

7:51:15 

5 

Yaw Data Gap 


706 

9;?,1:10 

9:32:05 

3 

Yaw Data Gap 


Subsequent 



3 

Yaw Data Gap 

230 

743-745 

0:00:00 

3:00:00 

3 

Normal (Mode 5) 


745 

2:38:32 



Start RCS 


745 

2:40:45 

2:59:10 

3 

Yaw Data Gap 


746 

4:21:00 

4:39:40 

5 

Yaw Data Gap 


747 

5:44:25 

5:50:10 

5 

Yaw Data Gap 


747 

6:01:30 

6:20:25 

5 

Yaw Data Gap 


748 

7:24:30 

7:31:10 

5 

Yaw Data Gap 


748 

7:42:00 

8:01:10 

5 

Yaw Data Gap 


748 

7:46:58 

7:48:22 


Burn 


749 

9:10:00 



Stop RCS 


749 

9:05:20 

9:11:30 

5 

Yaw Data Gap 


749 

9:22:40 

9:41:40 

3 

Yaw Data Gap 


Subsequent 



3 
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I'ablo. 4-4. Data Base Applicability for Maneuvers (Continuation 1) 


Event Times 


YAWINT Event 

Data Base Description 


1 to 819 

0:00:00 

6:10:00 

819 

6:17:50 

7:26:05 

819 

6:39:27 


820 

7:56:45 

9:05:55 

820 

9:20:36 

9:21:36 

821 

9:27:15 

10:49:45 

821 

10:12:02 


822 

11:18:10 

12:21:20 


Subsequent 


Subsequent 

Dp to 1072 
1072 

1071/1072 


0:00:00 24:00:00 

23:25:32 

22:32:20 23:54:10 



Normal (Mode 5) 
Yaw Data Gap 
RCS Start 
Yaw Data Gap 
Burn 

Yaw Data Gap 
RCS Stop 
Yaw Data Gap 


Up to 862 

0:00:00 

7:53:20 

4 

Normal (Mode 5) 

862/863 

7:41:55 

7:52:40 

4 

Yaw Data Gap 

863 

8:01:40 

9:00:25 

5 

Y’aw Data Gap 

863 

8:10:27 



Start RCS 

863 

9:22:22 

9:29:21 


Burn 

864 

10:00:02 



Stop RCS 

863/864 

9:22:40 

9:32:50 

5 

Yaw Data Gap 

864 

9:40:35 

10:39:30 

4 

Yaw Data Gap 


Normal (Mode 5) 
Start RCS 
Yaxi7 Data Gap 


Tahiti 4-4. Data Base Applicability for Maneuvers (Continuation 2) 


DOY 

Rev 

Event 

Times 

YAWINT 
Data Base 

Event 

Description 

From 

To 

233 

1073 

0:13{25 

1:38:55 

5 

Yaw Data Gap 


1073 

1:10:22 

1:10: 50 


Burn 


1073 

1:53:02 



Stop RCS 


1073/1074 

1:53:10 

3:15:50 

4 

Yaw Data Gap 


Subsequent 



4 
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SECTION V 


CHARACTERIZATION OF AD ACCURACY 


A. SUN SENSOR ERROR SOURCES 

Attitude determination in all three axes was affected by sun sensor errors. 
Pitch and roll were affected through the biases determined by c »mparison of Scan- 
wheel output data with concurrent sun sensor data. The achieved accuracies of 
bias determination are demonstrated in Section 3.2,2 of Reference 2-9. Yaw AD 
accuracies will be explained in the following paragraphs. 

The CSFC method of determining yaw attitude from sun sensor data is subject 
to three classes of errors: (1) coupling of roll and pitch AD errors into yaw 

AD, (2) sun sensor alignment errors, and (3) sun sensor observation errors. The 
error source mapping into yaw can be expressed as the first-ordei Taylor’s expan- 
sion of the yaw estimate: 


n 


Y 

e 


Y -f 
o 



i = 1 


9Y 

3(f). 

1 


Act). 


where 

Y = yaw estimate 
e 

Y = true yaw 
o 

6, = ith error source 

For this case, there are seven error source groups: pitch and roll AD error, 

misalignment of the sun sensor boresight in spacecraft clock, cone, and twist, 
and sun sensor observation errors in the two sun sensor orthogonal coordinates 
cone and cross-cone. The sensitivities (first partial derivatives) of yaw to 
these error source groups are derived in References 2-2 and 3-3 and are summar- 
ized here: 

d Yaw cos a 

_ s_ 

3 Rol!l tan 3 

s 


3 Yaw sin a 

s 

3 Pitch ~ tan 3^ 
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Yaw 


() «, 


= -1 


Yaw 


a (i 


sin (a^ - 
tan 3 


Yaw 




cos (x^ “ Ug) sin Bg 
tan 3 


- cos 3 


B 


Yaw 


S2' - 1 


J a 


ss sin 3 r^os 3„ cos 3 cos (a - a_) + sin 3 sin 3,.l 
si-B s sB s BJ 


Yaw 


3 3 

ss 


sin (a - a„) 
s B 


1 2 

tan 3 S« s B 

s J 


= cos 3n sin 3 cos (a - a„) - cos 3 sin 3„ 
2 B s s B s B 


S. 


= sin 3„ sin 3 cos (a - a„) + cos 3 cos 3r, 
B s s B s B 


where 


(t , 3 = clock, cone angles of sun in spacecraft coordinates 

s s 

fji , 3 , y clock, cone, twist alignment angles of sun sensor boresight 
in spacecraft coordinates 


(i , 3 = cross-cone, cone angles of sun in sun sensor coordinates 

ss ss 


B. AD ERROR BUDGETS 

References 5-1 and 5-2 give error budgets for Seasat attitude determination 
in all three axes. Reference 5-1 reflected the consensus of estimates from LMSC, 
GSFC, and JPL based on pre-launch expectations. Reference 5-2 was written after 
6 weeks of orbital operations, and reflected changes of some estimates based on 
this operational knowledge. Further operational experience and thought by both 
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JPI. ami CnSl'il luivo produmnl lurthor rol inonunU ol iho AD orror bvulyvHs* Tho (?SFC 
etU imatos iiiul tlu'ir rat louah* aro proiumlod in Soot Ion \ oi RoloVinua' 1?“9. Tho 
dPI. oat imatoa aro aumnuiri/.od in TaMoa and 

raMo ‘v 1, whioli t roat a o\ pit oh and ri>ll Ah aoouraoit'a, ha.s i\o\ow\\ 
d i 1 1 oroiii'oi; 1 riuii a similar t.dUi' in Kv'Um oium' Tlu' »illvuMti\>n loi' SoanwluH'l 

iioiao attor runninj\ aviMaj\i‘ liltorinp, had boon inoroaiunl io O.tbt iloj\ in KiM or- 
iMU'o but lurlbor oporat ional I'xporionoi' indiiMtod wood jn'i’tornunoo with tlu' 

IdO’.s piu iinl llltor, r.o t lu' allooation waa roturnod to (ho l^.Obdop, lovol, dSFi' 
analvaoa (Soot ion ol RtMoronoi' ahovHsI no Soanwhool .spood-dopondont 

iMUs'ls jMoator than tho 0,0d-di'}\ aoouraoy ol t ho i r analvtio motluul, ao that 
allooation waa nsluood to O.Od iloj\. Tho va Uun; lor ti'mporaturo variations and 
ihroahold atabilitv liad hi on inoroa.sod bv in RiMoronoi' ‘v-d to allow U^r tho 
poa^;iblv larju'r otlov'ta o\ oporat inj\ on a ainy,h' Soanwhoi'l ralhor than tho tlual 
UnniiinaM .swali'in; opi'ral ional oxpi'rionoi' has ^lhown no roaaon to airspool larp,or 
otloots, av> tho albn'alion^i woro ri'turnod to tho valuoa ol Koloronoo ‘v-1 . By 
similar riMaoniny, tho allooativ^n lor attitiulo oompul oi* drill wa.s rodnoovl lo tho 
\o\< ond ol tho ranpo o.ivon in Rolonmoo ‘■>--d . I’ho bias di‘i orminat ion orian* alio- 
v‘at ion waa roduood to tho valuo.s I'alimatod bv tlSFt' tor tho ilot orminat iona t Ivov 
porloimod. rho allooation tor random horiron radianoo oriaua** luui boon mado a 
lunotion ol aubsjxiooo i a I t lat itvulo, baaod on tlSl't' analwsia* tSootion l,a.d ol 
Koloronoo d Fianro ‘v^l ahow^^ t lu» t inio linu'tiiu'i lor tlu.s oi ror budvol and 

lor tho total Ah aoonrai*v ol pit oh am.i roll. 

It aluniUl bo iu>( od that Tablo S-d applio,s imly uiulor tho oonditioiu; whli'h 
booamo normal tor Soa,sal : ainalo Odaht'l Soanwhool oporation, with no aun intor- 

loronoo. rhoro woro poriodj? ol l imo durina dava ddO-^ddb t,aiul mnoh lator on 
rov 1 1 I wlu'u tho loll Soanwhool waa uaod; thoro waa an approximalo 0,S“doa 
oll.sot olosorvtiblo whon Soanwhoi'la woro ^iwitohod, and intormiltont .sun intorlor- 
onoo oauaod d-doo, orrora. On tho riaht Soanwhool I rom dav ddS to dav d‘vS, 
thoro waa an i\it or I iMonoo otloot onoo a rovolution Uwioi' a rovolulion 1 rom 
vlav dsn to dav dSS'), whii'h waa laraolv but not oomplotolv otlitod bv Iho tlSFt' 
data piiH'oa.sina j^y.stom. Alt hour, h *hta olloot voaa small aiul probably noal laH'^l^' 

it it happoiu'd durino^ Modo S, (h*ro woro aomo rovolutiona, amona I luMii rova /Oti-' 
,’11 plus 81d atul Si^l, whon Modv' s wa.s not uaod, and tho ol toots on t imtl AD wort' 
on I lu‘ tudi'r oi tl.d--l.ri 

Tablt' S-d j'ivoa t lu' ortan* Inulaot a \ov vaw dot t'rminat itui 1 rtuu .sun at'natu* 
data, Tlu' Inulaot a iov imlivldual t'rrtu* alnlrot'^^ tlillt'r 1 ivun t lunu' iM'ovituialv 
alltu'ait'd in Rt'U'it'uot' S-d in ihrt't' aroaa . Slnoo tomporaturo ai'ntlii'nta woro 
nunii'rali' otv.;ipa isul to t lu' maxima prt'tl 1 1 ' I t'd boltu'o lauiu'h, tlu' a 1 1 tu'a t 1 tui lor thoat' 
olli't'la waa rtuhu't'tl bv .1 laottu’ ol aj^piu^x ima t o 1 v d arul applii'd tmlv to t lu' axis 
about whit'h t lu' maximum t'l It'ot oould In' oxpootod io ooour. Al iannu'ut unot'iaalnty 
wa^. iiuM t'aai'tl 1 atnunint lor t lu' aoouraov ol Int liaht oalibratitm oi at'natn' 
lu'atl d am! I tu' t hi' appari'ut Innnula oi ml I 1 anmont ol aonatM* lu'atl 1 (aot' Si'o- 
t ion l.l.d oi Rt'loroiu't' d- at\d Subaoot itui 1 11 -lO . Tho oontributitm dno to pltoh 
ami It'll AD I'rrt'r waa nu'dilit'd 1 1 ' t’orroapoml It' aotual St'aaa. I'xpi'r I t'lnn' I'u I'at'h 
st'ust'r lu'atl. Tlu' It'tal vaw AD orrt'r va I uos alu'wn in d'ablo S-d art' t lu' rt'suit ol 
mappinp, t lu'^u' t'rrt'r >u'uitn'a intt' vaw with tho otjuat it'ua t'l Subaot'tit'n V-A. Tlu' 
ranao t'l valiu'ii at von i.s t lu' minimum to maximum valuoa t'bat'ivt'tl in aiu'h mappina 
ovt'i I lu' t'nlirt' ranat' t'l ram pt'.sitiona t'xpt'r i t'uootl bv Soaaat . 



Table S-1. AD Krror lhu4»i't s ( kO loi' Pit oh and Roll 
(lor t^poratioii on Ri^^ht Soamvhoel Only) 


Ni> 

. 

1 1 em 

Pit oil 

Roll 

1 


Random horir.on radianoe variation 

0.0S-0.P‘i 

O.OS-0. K’S 

> 


Ni'iiU' .11 tor U'O-H avorap.inp, 

0.01 

0.01 

1 

Soanwheel speed variat ion 

o.o:’ 

O.OP 

4 


TemjUM'at ure var lat ions 

0.0 (PI 

0.0 (A) 

S 


Thre.sho 1 d st ah i 1 I tv 

0.0 (P> 

0 .0(p) 

0 

Attitude oomputer drilt and ap.tnp. 

0 .00(1 

O.OtiP 

/ 


Bias determinat ion error 

O.Oh 

0.07 

8 

Hias dt‘tonuinat ion ohservahi 1 i ty 

0.01 

0.01 



TiU a 1 (dep, , lo) 

0 . 1 .'.-0 . 

0 . 1 1-0 . 1 8 



Tahle AD Krror lUuipet 

(5o) lor Yaw 




(When Sun is Observable) 


No. 


Item Sensoi' Head 1 

Sensor Head 11 

Sensor Head 5 

1 

Sun 

sensor aooiiraov 0.05 

o.os 

0.03 


(.pi'i- 

.ixiiO 



> 

Launoh shook 0.01 

0.01 

0.01 


Sun 

senso r all pnment 0 . 00 ) 

0.08'’ 

0.P3 


(.piM- 

ax is) 




Thermal ilistortion 0.0P5 (oone) 

0 . 0.15 U wi St ) 

0.0P3 (twi.4t) 

S 

IMto 

h and roll AD error 0.000-0. 

0.000-0. l.'i 8 

0.000-0. 1 OP 


ooupiinp, into yaw AD error 





Yaw AD error (dep, lo) 0,051-0.1105 

o.ooh-o.;i08” 

0.23 3-0. 108 


'\\ssnmos t hat known misal ip.mnont is taken into aoeount 
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ATTITUDE DETERMir4AT!Cr4 ACCURACY, deg (3a j 






Noto that Table ^)-2 accuracies for sensor head 2 apply only if correction 
has been made for the measured misalignment of this head. Such correction was 
made in GSFC data processiiag up to day 255. Definitive attitude files for days 
256“28‘3 can be corrected by using the information in Figure If the indicated 

correction is not applied to this data by the user, the accuracy estimate should 
be appropriately revised. 

The yaw AD accuracy fimction of time is presented in detail in Figures 5-9 
through 5“12, 

c:. YAW INTERPOLATION ACCURACY 

The accuracy of the yaw interpolation process can only be inferred, not 
measured. By definition, it is a means of supplying an estimate at times no 
direct measurement is available, so there is no direct means of verifying the. 
estimates made. Hc^wever, an estimate of interpolation performance can be formed 
by comparing its performance against real data under simulated data gap conditions 
and making appropriate inferences. 


1. Data (Jap Simulations 

a, Qu 1 e s c en t Per iod . Data from three consecutive days (214-216) of the 
first quiescent period, 37 revolutions in all, was used to simulate the attitude 
data conditions of the second quiescent period (days 256-283). The parameters of 
YAWINT data base 2 (see Table 4-2) were used to create an interpolated yaw esti- 
mate funct ion for the orbital interval between -1040 and +4100 s from ascending 
node. This Interval approximately covers the orbital phases of sun sensor data 
gaps in the second quiescent period (see Figure 2-12). To the extent that yaw/ 
roll behavior is similar in the two quiescent periods, this simulation models 
yaw interpo lat ion for days 256-283. Three consecutive days were chosen to aver- 
age out the effects of longitudinal variation, since Seasat was in a near-3-day- 
repeat orbit. Figure 5-2 presents the results of this simulation as mean error 
(lO, mean plus and minus standard deviation of error (p ^ iO , and the extrema of 
error, all as functions of time from data gap start. Note that true yaw tends 
to be more positive than the interpolation in the northern hemisphere and more 
negative than the interpolation in the southern hemisphere. The average mean 
error over the entire simulated gap is -0.01 deg. Without corroborating evidence 
f rom an independent source (e.g., the SAR) of data from the day 256-283 period, 
the dynamic function of error has low confidence as a predictor for that period, 
and the mean error slunfLd be considered zero. Overall performance is 0.79 deg 
(3 RMS). 

To estimate the similarity of actual (observable) yaw in the target era 
with simulated interpolated yaw, another comparison was made. The RMS error of 
interpolated yaw in the orbital phase 4100 to 5000 s from ascending node was 
evaluated for d/iys 214-216 in the first quiescent period and for days 257, 258, 
and 268 in the second. Sun sensor data does generally exist ‘ in this orbital 
phase for the second period, and, therefore, a real performance estimate can be 
made. As it happens, this orbital phase is quieter (in yaw) than that 
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DATA GAPS SIMULATHD FOR 
CECCrJD QUIESCENT PERIOD 
YA// & ROLL DATA FROM DAYS 214-216 
YA//ir;T DATA BASE 2 



TIWE FROM DATA GAP START^ 





ri'privuMit ihI In lirinr ' - at h\uU lor vla\^i .M f- iho oompari.^u>n o\ ro^ullt^^ 

botwoiMi t ho t il^;l aiu! ^u^^Mul viuiosi‘ont porituis i ^ihown Kitoi in labU' n 
riuMH' ia an inoroaso in RMS oi’ror lor lliis orbital pirnn' i rom t ho l irnt poriiul 
to t ho jJOi-i'iul. iioatinp. thi^‘> inoroano as duo io a provivui.alv iinav oount od indo~ 
pondont orriU' ^iol^Ol’, it woiiKl havo a valuo ol tK > ' dor, RMS'). I’lii^s oan bo 
KSS^tui with I ho dop, O RMS'i ovorall portonnanoo in the' r.imulatod data paps 

to inlor dov» t, ^ RMS) lor t ho roa I data r,aps in t iio stH'oiul quiosiont porioil, 

b. M’^^'do a Tor i od . As dosoribod abovo lor t lu' tpiioscUMit poriods, data 
p,tips \v’oro s i mu I at od i A^r t lu' Modo poriods, Sinoo roai liata paps had alro*idy 
appoarod at t ho lu'pinninp, oi this tu'a, onlv liinilod data was availablo lor ooni- 
pari son. i’rom dav dd"* to day d hi , tho aotiuil data paps variod ovmks i dorab 1 y in 
nunibor, .siao, aiul phas.o, as k‘an bo si'on in l’ip,uro d-ld, Wy day d , thori' was 
onlv ono p.a.p, aiul it v;as botwoon ^HU) and aO‘U^ .r I rom asoiMulinp, nodm I’sinp, 
YAWINI’ data loiso \ and tho av.ulablo data I rom davs d!\S aiui dd^\ tho day d bV 
data p,ap was simulutod, I’ho ro.su 1 1 s aro slunvui in I'ipairo S- n Tho pal torn o[ 
orror is i^bvioirslv ditlA'ront 1 rom that ol tho quiosoont pi'iiod U'Mpuro b-d'l, 
drivon mainlv bv t lu* Miulo 'i oxoursion in t !u* lattor part ol t lu* pap. Tho ovoim 
moan ij; about O.iM dop and ap:iin should bo oonsidorod aoro. 

Uooairso oi tho s i pn i 1 i oant 1 v pooror i nt o r['o 1 a t i on por l oimauoo noar t lu* 

Modi* t'Xt’ursion i,whioh was in t lu* sout horn lunn i splu* lu* , lu*ttc*r i*stimatos ol 
pi*r lAu'maiuu* i*an In* madv* by oonsidorinp tho Modt* rt*pi(Ui as. distinot 1 rom tho 
rost oi t lu* rovt'lulion. Aovum\ 1 i up, 1 v , i n t orpo 1 a t i on pta* tAnamnuu* was ostimatod 
for vlata in tho intoinal of IdlR) s t'ollovdnp tho ri'll ^^ip,^aI disoonnoot (,Modo 
starl'l aiui si'paratolv tor the* data in t lu* rost ol t*aoh rovolution. Tablo A- ^ 
shows tho I'osults lor vaw prodiotion CVp"! from tho I'ouric'r Si*rii*s, yaw intorroil 
from roll and I i rati* O' 0 ^ ^^nd t lu* ovorall yaw i n t orpo 1 at i on oslimato O’^ . 

As Ovin bo soon in tho tabl,*, i n t orpo 1 a t i on jn* r i o rnuiuv'O is t wi> to throe* t imos 
bolt or away T rom tho Modo a intorval than iu*ai’ it . I'his mtoms that soiontilio 
data 1 rom tho iu*rt lua n homisphoro is muoh bot t or sorvod than that I'rom tho 
soutlu'rn homispluMo; iiappilv, tho lUM'thorn homispiioro a Uu> liraws tho most 
so i i*nt i t i o i nt oro.s t . 

l\*r I o rmaut'o ovaluations wore* mado t rom availablo data t hrouphout t lu* day 
dd‘>~ddS iH'riod. I'oi' oaoh individual day within that porioil, tho appropriato 
aoouraoios oan bo obtainovl bv i nt orpo I at i up Tablo TIu* data Trom days 

and dda tt*iuls to In* oonoonlratod in tho low latitudo^^, so that thoir inrlormatuu 
1 ipuros aro biasod toward hip,lu*r unoo rt a i n.t i os . A moro appropriato st*t of 
valuos lor i lu* davs lu*twoon .MO and d O'* aro l lu* ovorall i ipairos p,ive*n by tho 
tiiblo tor this poriod. 

Tho orbital pliast* oT tho data paps in tho lator portion of tho Mode* A 
pi*riod \v’as alnuust tho .‘oimt* a^^ Tor tho st't'ond ijuiosoont iH*riivl. Suoh data piips 
Wi*ro simulati'd, ap.ain on data 1 rom davs and but us i up, YAWINT datti 

baso s. The* rosulls aro shown in IMpuiu* , Tlu* ovorall moan orror is now 

O.thS'dop,, whioh mav bo sipnitioant. That tliis valuo i .s larpor than lor tlu* 
otlu*r poriods i duo primarilv to t lio pooror |H*rlormanoo noar tho Modo S oxour- 
sii>n. Tho i n t oriio 1 a t od ostimato oonsistoiUlv undorshot t ruo vaw ior larpo 
oxoursions, in ordor to minimiro orror Tor tho moro usual small oxoursions. Tlu 
si. a* o! tbit-; midorsh^'ot was as lar u* as d dor,. 



Tabk' Yaw Iiitorpolat ion Accuracy (Muivlie 


of a Long Data CJap) 


Ver fornuiiu'e (3 x I<MS, dtHD for Givon »rl^ital Phaso and YAWINT Data Base 2 


Data 


Sun Sensor Data 

Available 

Sun 

Sensor Data Not Available 

Day 


in Tarp,et 

Kra 


in Target Bra 



AY 

P 

AV., 

A 

AY 

AY 

P 

AY 3 

A 

AY 

214 

0 .t )8 

0.37 

0.58 

0.82 

0.95 

0.77 

213 

0.71 

0 . 3b 

0.385 

0.82 

0.95 

0,77 

21 b 

0.7 3 

0.53 

0.58 

0 . 86 

1.00 

0.83 

Overall 

0.71 

0.33 

0.58 

0.83 

0.97 

0.79 

237 

0.73 

0.78 

0 . 80 




238 

0.8b 

0.73 

0.73 




2b8 

1 . 04 

0.80 

0 . 82 




Overall 

0.89 

0.77 

0.79 

0.99‘^ 

1 .ir* 

0.95^* 




MoUi> ■) 

Start + 1200 

s 

NOT 

(Mode 5 

Start + 1200 

s) 

Data 

Day 

AY 

P 


A 

AY 

(d.b. n 

AY 

(d.b. A) 

AY 

P 

AY,, 

A 

AY 

(d.b. 3) 

AY 

(d ,b .4) 

22b 

1.81 

1.82 

1.47 


0 . 9 3 

1.02 

0.68 


228 

' l.b2 

1 . 94 

1.45 

1.79 

0 . 82 

0.91 

0.62 

0.81 

229 

1.79 

2.01 

1.67 

1.90 

0.8:^ 

0.93 

0.63 

cc 

0 

2 32 

2.4b 

1 . (>0 

1,85 

1 

0.83 

0.72 

0.63 


2 34 

3.42 

2 . 04 

2.5b 


1.14 

1 . 10 

0.97 


Overal 1 

2.32 

1.89 

1 .84 

1.84 

0.92 

0.94 

0.72 

0.82 

230 

3.84 

1 .9 1 


1 .89 





251 

4. 18 

2. 19 


2.22 





Overal 1 

4,01 

2.05 

NA 

2 .06 

3.7r^ 

1 . 19 '’ 

NA 

1.24“ 


‘hnferrod from comparison of avaihiblo data from appropriate orbital and mission phases. 


NOTES: 1 . DATA FROM DAYS 228-229 

2. YAWINT DATA BASE 3 







In the lator \nu'tion of t ho Modo ora, tho only yaw bohavior ohsorvahlo 
by tho sun sonsor was in tho orbital phaso dlsfurbod by Modo Comparison of 
avillablo data from that t Imo poriod with data from days 22^-229 havin ^5 tho samo 
phaso rolationship to Modo 5 start shows a dooroaso in intorpol at ion porlormanoo 
from K84 doy, to 2 *06 doy (3 RMS). Pol hawing tho rationalo usod fi^r tho qulos- 
oont poriod, this additlcnial unoortainty oan bo attrlbutod to anothor orror 
so ur 0 o o f va I uo 0 • 9 4 d o y ( 3 RMS ) . Comb in ing t h 1 s in an RMS s o ns o w i t h t ho 

0.82-dog t>vorall porformanoe for days 228-229 loads to tho ostimato of 1.24 dog 
(3 RMS) for yaw Intorpolat ion porformanoe outside tho Modo 5 disturbance region. 

2. Intorpolat ion Error Distribution 

Tho i ntorpo lat ion orror of Figures 5-2 through 5-4 i^as analyzed to dotor- 
mino tho relative frequency of various values of orror. Figures 5-5 and 5-h are 
the histograms resulting from this analysis, and incorporate 882 and 492 samples, 
respectively. The only truly significant difference between the two histograms 
is the mean shift due to Mode 5. The apparent deviations from tho classical 
bell-shape are probably due to tlie small number of samples. 

Figure 5-7 shows the distribution of RMS error with respect to orbital 
phase for the three simulated interpolation cases. As for the previous two 
figures, the scaling is correct only If true yaw behavior in the target (simu- 
lated) eras was identical to that in the reference eras. However, even if 
actual behavior is only similar to these simulations, it appears that RMS error 
is less in the oibltal phase between +72^^N and descending node than in the rest 
of the revolution. This information could be used to give preference to data 
taken in this orbital phise, but was not used in tlie overall estimates (given in 
Subsection V-C-3) of yaw interpolation performance. 

Figure 5-8 plots the cumulative distribution function of RMS interpolation 
w^rror for four cases: (1) simulated second quiescent period, (2) simulated 

early Mode 5 period, (3) simulated late Mode 5 period, and (4) a Ckuissian dis- 
tribution. Each data set was normalized to its ensemble RMS value to help the 
comparison of distribution shapes. It is quite evident in Figure 5-8 that 
there is no significant difference between ciny of the functions and the normal 
distribution. The raggedness at high multiples of RMS is due to the paucity of 
samples. Therefore, RMS interpolation error can be understood as the standcird 
deviatirn of a normal distribution, so that the ninety-ninth percentile is at 
2.56 RMS, the fiftieth percentile is at 0.67 RMS, etc. 

3. Yaw Interpolation Accuracy Versus Time 

Figures 5-9 through 5-12 reflect all of the information in the preceding 
paragraphs, and they display yaw attitude determination accuracy as functions 
of orbital and mission phases. Figure 5-9 shows the early mission phases when 
the sun was continuously observable by sun sensor head 1. The variations in 
yaw AD accuracy are caused by the pattern of sun positions within the sun sensor 
FOV on the given day, but without regard to attitude deviations. Significant 
yaw excursions did occur previous to day 195 and after day 219; the former 
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SIMULATED FROM DAY 214-216 
AND DAY 228-229 DATA 



SECOND QUIESCENT PERIOD 
early MODE 5 PERIOD 
-- later MODE 5 PERIOD 


S = RMS (YAW - INTERPOLATED YAW) 




PERCE 



M (MULTIPLE OF RMS YAW INTERPOLATION ERROR) 

Figure 5-8. Significance of Yaw Interpolation 
Error RMS Value 
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would cause no discernible change In Figure 5-9, but the latter would cause a 
shift fon«ird in apparent mission phase of approximately a day for each positive 
degree of yaw. There were also some 15- to 20-min periods in several revolution^ 
betwecti days 220 and 225 when pitch and roll were seriously in error due tii sun 
interference; the coupling of pitch and roll Al) eri'or into yaw could produce 
1 to 1.5 deg of yaw AD error at those times. Allowance for these variations in 
accuracy would have to be done on a rev-hy-rev basis, which is impractical in 
this report. 

Figure 5-10 covers the mission period in which yaw data gaps began to 
appear and the transition to sun sensor head 2 occurred. The interpolation 
accuracy as a function of orbital phase was computed with the variance model 
described in I’ubsection Il-B-3, using variance values inferred from Table 5-3. 

One variation from that model is the use of the Table 5-3 value for AY perfor- 
mance as the limiting variance value rather than the less conservative oyo com- 
puted by the model. Again, actual attitude deviations were not taken into 
account in this figure, so that the true pattern of accuracies can differ from 
that shown; l\owever, Figure 5-10 is still a good guide to the accuracies avail- 
able in this mission pliase. These accuracies do not apply to those revolutions 
on days 227, 230, 235, and 238 when the RCS/gyro control mode was in force. 

Except for the half-rev transition (each way) between control modes, tlie yaw 
interpolat ion accuracy in the RCS/gyro mode was 0.7 deg (3 RMS). 

Figure 5-11 shows the later portion of the Mode 5 mission period. All 
directly observed yaw data comes from sun sensor head 2. Accuracies were com- 
puted as described earlier for Figures 5-9 and 5-10, except that an additional 
uncertainty was inserted in the orbital phase during which Mode 5 was used in the 
northern hemisphere (see Subsection lII-C-4). 

Figure 5-12 covers the second quiescent period and includes the end of the 
mission. The .scale is expanded by a factor of two compared to Figures 5-10 and 
5-11 because of the smaller range of uncertainties to be expressed. Note that 
sun sensor head 3 produced data over the high northern latitudes during the last 
few days of the mission, starting around day 274. 
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SECTION VI 


SUMMARY AND RECOMMENDATIONS* 


A. SUMMARY 

The first satellite built to observe the world’s oceans, Seasat suddenly 
stopped sending data from orbit on 10 October 1978 after 105 days in space. In 
October, it was believed that a massive short circuit of the solar array 
rendered it unable to transmit or receive any signals. It was concluded finally 
that the loss of power was caused by a massive and progressive short in one of the 
slip ring assemblies. (The slip ring assemblies were used to connect the power 
system in the Agena with Seasat ’s rotating solar cell panels.) 

Although the premature termination of the Seasat mission prevented observa- 
tion of the change in sea conditions from season to season, a sufficient amount 
of flight data was obtained to evaluate the feasibility of the scientific aspects 
of the mission and to perform comprehensive post-launch attitude determination 
evaluations of the following: the Earth’s horizon radiance model, cold cloud 

effects on IR data, the performance and the alignment accuracies of the attitude 
sensors and control system, and the capabilities of the Seasat attitude support 
software. Furthermore, the results of these analyses suggested the following 
recommendations for future mission support. 


B. RECOMMENDATIONS 

The following recommendations for attitude determi ation mission support 
hardware and software are derived from the Seasat mission support experience. 

The hardware recommendations deal primarily with the IR scanners and sun sensors, 
and some additional attitude sensing hardware is suggested. The software recom- 
mendations deal briefly with each aspect of the data processing to point out some 
of the handicaps in the operational support software. No comments are presented 
about the observed sun interference in the IR scanner. 


1. IR Scanner Recommendations 

The attitude determination error analysis demonstrated that the increased 
accuracy that was expected from the normalized threshold locator logic was nulli- 
fied by the relatively broad band pass of the IR scanner optics. As a result of 
the broad band pass, atmospheric phenomena on the Earth’s surface as low as 8 km 
(4.3 nm) Ititude perturbed the threshold normalization voltage. The random 
errors thus introduced could not be modeled and exceeded the horizon radiance 
corrections being applied by the attitude determination software by a factor of 

2, Although the normalized threshold system may have been necessary for attitude 
control stability (in that it reduced long-term fluctuations due to the 
seasonal/latitude-dependent radiance variations), a fixed threshold system, 


*Based heavily on Section 5 of Reference 2-9. 
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assuming the same IR band pass, would have presented a more tractable problem 
for attitude determination. For the fixed threshold system, unmodeled cold 
cloud errors would have been greatly reduced and the errors in the model for the 
systematic radiance changes would have been no worse than those for the normalize 
threshold locator. Therefore, the attitude determination errors for the fixed 
threshold locator logic would have been smaller than those for the normalized 
threshold system. 

The attitude error analysis completed 1 year before launch demonstrated 
that the 3o error limits on attitude determination performance could not be met. 
The analysis was initiated and completed too late in the mission support to have 
any effect on the spacecraft flight hardware configuration. Although similar 
analysis performed by the spacecraft contractor much earlier in mission support 
had concluded that the error budget could be met, it is recommended here that 
such early analysis be performed by a third party or by the party responsible for 
attitude determination data processing. 


2. Sun Sensor Recommendation 

Ciiven the problems that occurred with sunlight in the IR scanners and that 
resulted in control anomalies that reduced the effectiveness of the yaw inter- 
polation algorithm, more comprehensive sun data coverage should have been provided. 
The six-sun sensor configuration proposed by GSFC before launch would have been 
invaluable for the determination of yaw, given the control problems that existed. 
Although sun sensor coverage near the zenith was discouraged because unfavorable 
geometry caused degraded accuracy, the addition of two or more sun sensors in 
this area would have provided increased sun data coverage near zenith, shorter 
sun data gaps, and a higher overall determination accuracy. 

During the spacecraft design review, a request was made for information 
about the behavior of the sun sensors at the times when the sun traversed the 
edges of the FOVs. If this information had been communicated, significant pro- 
blems in the data processing at these times could have been avoided. Software 
specifications to specifically eliminate or analytically correct anomalous data 
from the sensors, near and beyond the edges of the FOV, could have been implemented 
The obvious advantages would have been efficient and uninterrupted attitude data 
processing and possible extension, with some degraded accuracy, of the effective 
sun sensor FOV using the fine reticle telemetry data to angles as high as 60 
deg from the boresight. 


3. Yaw Interpolation Hardware Recommendations 

Another alternative to the problem of yaw attitude determination and yaw 
interpolation involves the use of the yaw gyro as an attitude determination device. 
The yaw gyro was a critical instrument for yaw attitude control using hydrazine 
thrusters during orbit adjust maneuvers. For this reason, in addition to the 
fact that the gyro had a specified operating lifetime of no more than 30 days of 
continuous operation, the use of the gyro data in the attitude determination 
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system was negated. However, given the problems with yaw interpolation caused 
by the sun- interference-induced control anomalies, periodic use of the gyro data 
would have been invaluable for testing the yaw interpolation parameters and for 
evaluation of the interpolation absolute accuracy. In the pre-launch phase, 
during the attitude software specification process, some consideration was given 
to include the gyro telemetry data in the attitude telemetry record, but this 
data was excluded from the record because of space limitations on the input disk 
data set, and because there was low probability that the gyro would be available 
for this purpose. Magnetometer data, however, were included, even though the 
relative accuracy of this data, limited primarily by the resolution of the tele- 
metry word, was known to be insufficient for mission support. 

It is, therefore, concluded that during early mission support, a method of 
checking the yaw interpolation parameters at all phases of the mission should 
have been considered. A schedule of weekly or twice-weekly one-orbit runs with 
the yaw gyro "on,” during low sun data coverage, could have been a nominal mission 
support procedure. A second backup yaw gyro package may have been required for 
this support, and, in the event that control anomalies had occurred, more frequent 
use of the gyro could have been justified. 


4. Magnetometer Recommendations 

Following the discovery that the yaw interpolation accuracy was being 
significantly degraded by the spacecraft control anomalies, magnetometer data 
was considered as a possible source of more accurate yaw attitude information. 
However, the following problems were encountered: 

(1) The resolution of the telemetry word (4.4 milligauss) was insufficient 
and restrictive for accurate yaw attitude, particularly at the high 
latitudes when the horizontal field was low. 

(2) The magnetometer triad was ground-calibrated against the sun data, on 
the full sun orbits, using a least-squares method, and unusually 
large biases were discovered along with evidence for sizable misalign- 
ment, crosstalk, or nonorthogonality of the triad. 

It is, therefore, concluded that the telemetry word size for the magnetometer data 
could have been increased to be comparable to the noise in the magnetometers. A 
method of electronically nulling the biases on each axis of the magnetometer could 
have assisted the control system in dumping angular momentum. More accurate pre- 
launch calibration data from the manufacturer, or pre-launch calibration of the 
device once mounted on the spacecraft, would have been desirable. 


5. Software Recommendations 

Operationally, the Seasat Attitude Determination System (ADS) performed well 
after it was corrected for software errors and operational inconvenience. The 
Seasat ADS was composed of the following four subsystems: (1) Telemetry Processor 

(TP), (2) Definitive Attitude Determination System (DADS), (3) Log Interrogation 
and Data Management Utility (LIDM) , and (4) Yaw Interpolation Utility (YAWINT) . 
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The software was specified to be modular, with separate load modules for 

each major data operation. In view of the Seasat experience, this remains a 

recommended design feature if the number of subsystems required for daily mission 
support is small. S. zsat operations personnel were able to learn the four sub- 
systems required for attitude support, and analysis personnel were able to write 
detailed operations methods for four subsystems without difficulty. 

The DADS could have been used only to apply corrections for biases, oblate- 
ness, and horizon altitudes. If no smoothing option had been applied to the data 

in the DADS subsystem and ephemeris-related operations could have been performed 
only at the output period, a consider ible increase in the speed and volume of dati 
processed in a single pass through the DADS could have been realized. Alternative 
the array sizes required for data processing could have been reduced, possibly 
resulting in a maximum core requirement of 400K bytes for 30 min of data. 

YAWINT, which performed yaw interpolation, could have been expanded into an 
output data smoother, incorporating v/hatever smoothing and interpolation method 
was necessary for all three functions (pitch, roll, and yaw). Methods such as 
Chebyshev fitting or spline fitting procedures to smooth the data and bridge the 
gaps in pitch and roll could have been specified as options. However, after view- 
ing the repetitious pitch and roll functions from orbit to orbit throughout a 
given day, it became apparent that another method would have worked well; the 
average pitch and roll versus mean anomaly, for a given day’s data, could have 
reliably bridged gaps of up to one orbit in length. 

Computations for the log (such as minimum and maximum values) could also 
have performed by this subsystem. The logging functions specified for Seasat were 
extensive and comprehensive in their goals. During the early post-launch phase, 
no daily routine was established for data processing mission support, and because 
of excessive problems with the overall mission, analytical and software personnel 
were unable to monitor and quality assure the logging aspects of data processing. 
Therefore, the log was neglected and never met (nor was required to meet) the 
expectations of the specification. If the mission had continued, with time 
allocated for log functions maintenance in all subsystems, the log would have 
been invaluable in monitoring and perfecting the data processing procedures. 
However, for this mission, as for most missions in the early orbit phase, only a 
printed page of statistics for each day processed is necessary for the log. If 
this option had been available and no operational interface between the log and 
the operator would have been necessary, early mission procedures could have been 
simplified. Version 1 of the Seasat software could ve been less complicated, and 
acceptance testing could have been simplified. Also, mission support programmers 
working on a later version of the software (scheduled for completion 3 to 6 months 
after launch) would have been ready and available to assist with software-related 
processing problems. 

Further specific ground software recommendations are provided by GSFC in 
Section 5 of the Reference 2-9. 
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APPENDIX 


AD 

ADS 

AOT 

CLA 

DAD 

DADS 

DAF 

DOF 

EOM 

FOV 

GSFC 

IPD 

IR 

JPL 

LIDM 

LMSC 

OACS 

PDFS 

QA 

RCS 

RMS 

RSS 


ABBREVIATIONS MID ACRONYMS 

Attitude Determination 
Attitude Determination System 
Attitude Orbit Tracking (Tape) 

Control Logic Assembly 

Definitive Attitude Determination 

Definitive Attitude Determination System 

Definitive Attitude File 

Definitive Orbit File; Degrees of Freedom 

End of Mission 

Field of View 

Goddard Space Flight Center 
Information Processing Division (GSFC) 

Infrared 

Jet Propulsion Laboratory 

Log Interrogation and Data Managei^ent Utility 

Lockheed Missile and Space Company, Inc., Sunnyvale, CA 

Orbital Attitude Control System 

Project Data Processing System 

Quality Assurance 

Reaction Control System (LMSC) 

Root-Mean-Square 

Root-Sum-Square 
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SAR 

Synthetic Aperture Radar 

SARA 

Satellite Alignment Reference Axe 

SASS 

Seasat Scatterometer System 

SSM 

Sensor Support Module 

TMDF 

Telemetry Master Data Flit 

TP 

Telemetry Processor 

w. r . t . 

with reference, to 


YAW I NT 


Yaw Interpolation Utility 


